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1  INTRODUCTION 


At  the  outset  of  this  work,  a  controversy  had  begun  concerning  the  near-millimeter-wave 
(NMMW)  transmission  of  the  atmosphere.  Observations  of  attenuation,  then  and  since,  have 
been  reported  [e.g.  1-3]  which  are  higher  than  anything  attributable  to  accepted  theory. 
Absorption  in  this  region,  between  90  and  1000  GHz,  comes  primarily  from  rotational  lines  of 
water  vapor.  While  it  is  difficult  to  predict  absorption  in  the  wings  of  lines,  the  referenced 
observations  were  particularly  anomalous  in  one  or  more  respects: 

1.  Attenuation  per  unit  water  vapor  concentration  seemed  to  increase  with  decreasing 
temperature  much  faster  than  expected.  Energy  levels  in  the  water  molecule  and  expected 
intermolecular  interaction  energies  could  not  account  for  the  results. 

2.  Absolute  attenuation  coefficients  were  very  large  at  frequencies  between  water  lines  where 
low  attenuation  was  expected.  Instances  of  this  occurred  during  conditions  of  high 
humidity,  fog,  and/or  low  temperatures.  Values  of  extra  absorption  around  10  dB/km  or 
higher,  if  correct,  would  have  an  adverse  effect  on  applications  of  the  NMMW  band. 

3.  Discrete  absorption  features  appeared  which  were  not  at  frequencies  of  known  water  lines 
or  of  other  known  atmospheric  constituents.  Spectroscopic  studies  of  pure  water  vapor 
under  carefully  controlled  laboratory  conditions  have  not  shown  such  features,  thereby 
eliminating  possible  equilibrium  concentrations  of  water  dimers  as  a  potential  explanation 
[4,5]. 

Although  oxygen,  ozone,  and  various  trace  species  also  absorb  in  this  region,  such  phenomena 
cannot  be  assigned  to  them  and  have  been  termed  "anomalous  absorption".  It  has  been  suggested 
that  exotic  polymers  of  water  might  be  responsible  [2].  On  the  other  hand,  contrary  results  have 
also  been  reported  [e.g.  6,7],  and  assertions  made  that  there  were  errors  in  those  observations  in 
which  anomalous  absorption  appeared.  These  matters  are  covered  in  greater  detail  in  a  review  [8] 
published  during  the  course  of  this  work  and  included  in  Appendix  D  of  this  report. 

The  purpose  of  our  current  work  was  to  prepare  new  instrumentation  for  improved 
observations  of  atmospheric  absorption  under  conditions  purported  to  be  associated  with 
anomalous  phenomena.  We  chose  a  spectroscopic  technique,  rather  than  single  frequency 
transmissions,  so  we  could  look  for  anomalous  absorption  features  like  those  described  in 
above.  The  ability  of  the  atmosphere  to  change  significantly  during  time  intervals  of  the  order  of 
a  few  minutes  can  have  a  significant  effect  on  accuracy  in  recording  spectra.  Distortions  due  to 
atmospheric  changes  were  one  of  the  possible  sources  of  errors  cited  by  critics  of  anomalous 
absorption.  This  consideration  was  a  primary  factor  in  our  instrumental  design. 

There  are  two  basic  approaches  to  spectroscopy  in  the  NMMW  band.  Spectra  can  be 
obtained  from  tunable  narrow-band  sources.  Due  to  their  relatively  high  power,  good 
signal-to-noise  ratios  can  be  attained  at  any  given  frequency.  However,  there  is  often  a 
significant  amount  of  time  and  effort  devoted  to  tuning  so  that  it  is  difficult  to  get  a  high 
resolution  spectrum  in  a  time  which  is  short  compared  with  that  taken  by  some  changes  in  the 
atmosphere.  The  alternative  is  Fourier  transform  spectroscopy.  The  disadvantage  here  is  that 
wide-band,  black-body  sources  are  relatively  weak.  Even  with  a  sensitive  liquid-helium-cooled 
detector,  long  integration  times  are  required.  If  a  long  time  is  spent  recording  each  sequential 
point  of  the  interferogram,  then  the  total  record  time  may  be  minutes  to  hours.  During  this  time 
any  changes  in  atmospheric  transmission  will  have  the  effect  of  a  modulation  on  the 
interferogram  and  will,  to  some  degree,  distort  the  spectrum,  which  is  obtained  by  Fourier 
transformation.  However,  all  frequencies  are  received  at  all  points  in  the  interferograms  recorded 
tthe  multiplex  advantage),  and  rapid-scanning  techniques  to  obtain  interferograms  in  a  few 


seconds  have  been  devised  for  use  at  infrared  wavelengths  [e.g.  9,10].  When  integration  times  of 
a  few  minutes  are  required  to  build  up  the  signal-to-noise  ratio,  multiple  interferograms  or 
spectra  can  be  averaged.  Then,  although  the  result  is  an  average  over  that  period  of  time,  each 
spectrum  in  the  average  is  practically  undistorted  by  atmospheric  change. 

We  elected  to  extend  rapid-scanning  FTS  techniques  to  the  NMMW  band  as  the  best  way 
to  look  for  all  three  purported  aspects  of  anomalous  absorption.  In  addition,  we  have 
incorporated  polarizing  beam  splitters  in  the  Fourier  transform  spectrometer  (FTS)  which  offer 
much  improved  efficiency  over  ones  used  previously  in  transmission  measurements.  The 
complete  system  described  in  this  report  is  designed  to  make  transmission  measurements  over  a 
novel  two-way,  400  m  path.  A  liquid-helium-cooled  InSb  detector  with  a  short  response  time  is 
incorporated.  Partial  tests  of  the  instrumentation  will  be  described,  but  difficulties  with  the 
detector  have  made  it  not  possible  to  obtain  atmospheric  transmission  measurements  as  of  the 
time  of  writing. 

2  FOURIER  TRANSFORM  SPECTROMETER 
2.1  Selection  of  Interferometer  Type 

A  two-beam  interferometer  (like  a  Michelson  interferometer)  is  an  efficient  way  to  obtain 
interferograms  from  which  spectra  of  signal  intensity  are  derived  by  Fourier  transformation. 
Unlike  a  grating  spectrometer  with  slits,  the  input  and  output  apertures  can  be  large  to  promote 
throughput  of  radiation. 

The  history  of  interferometer  development  in  the  NMMW  band  is  marked  by 
improvements  in  the  beam  splitter.  An  ideal  beam  splitter,  divides  the  input  beam  into  two  beams 
of  equal  intensity.  Thin  metal  films  cannot  be  used  here  due  to  absoiption.  Michelson 
interferometers  first  appeared  with  thin,  plastic  film  beam  splitters  [11].  However,  such  films 
had  less  than  ideal  reflectivity,  absorption  losses,  and  troublesome  gaps  in  efficiency  caused  by 
interference  within  the  films  [12].  Moreover,  the  angle  of  incidence  of  radiation  to  the  films  was 
often  close  to  Brewster’s  angle,  so  one  plane  of  polarization  was  not  effectively  split.  In  1969, 
Martin  and  Puplett  [13]  introduced  a  type  of  interferometer  that  also  restricted  throughput  to  one 
plane  of  polarization  but  incorporated  a  beam  splitter  with  nearly  ideal  performance. 

This  was  achieved  with  a  wire  grid  beam  splitter  that  reflects  one  plane  of  polarization  and 
transmits  the  orthogonal  one.  The  reflection  and  transmission  of  the  two  polarizations  approach 
ideal  if  the  period  of  the  grid  lines  g  is  much  less  than  the  wavelength  X.  (See  Appendix  B.) 
When  the  grid  lines  are  oriented  at  45°  with  respect  to  the  input  polarization,  such  a  beam  splitter 
can  divide  the  radiation  into  two  beams  of  nearly  equal  intensity.  While  early  grid  beam  splitters 
were  made  by  stretching  fine  wires  across  a  frame,  we  have  since  found  [Ref.  14  and  Appendices 
A-D]  that  high  performance  can  more  reliably  be  obtained  with  a  pattern  of  metal  strips  etched 
on  a  thin  plastic  substrate.  When  the  substrate  is  very  thin  with  respect  to  X,  it  is  essentially 
invisible.  Beam  splitters  of  this  type  are  used  in  the  system  described  here. 

Many  versions  of  the  Martin-Pupiett  polarizing  interferometer  have  been  described,  two  of 
which  are  illustrated  in  Fig.  1.  Part  a.  shows  a  version  with  separate  input  and  output  beams  like 
those  in  a  traditional  Michelson  interferometer.  The  output  beam  takes  a  different  direction  due 
to  the  rotation  of  polarization  caused  by  the  dihedral  mirrors  [13].  Part  b.  shows  what  happens  if 
the  mirrors  are  planar:  the  output  radiation  goes  back  in  the  direction  of  the  input  and  must  be 
separated  out  with  another  polarizing  beam  splitter  used  as  an  analyzer. 
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Fig.  1  Types  of  Martin-Puplett  polarizing  interferometer.  The  grid  polarizers  have 
their  grid  lines  at  45°  with  respect  to  the  page.  The  input  plane  of 
polarization  is  out  of  the  page,  a)  using  dihedral  mirrors,  b)  using  plane 
mirrors. 


While  separation  of  input  and  output  beams  is  a  nice  feature  of  the  layout  in  Fig.  la, 
dihedral  mirrors  are  difficult  to  align  optically.  One  axis  of  tilt  and  the  orientation  of  the  dihedral 
roof  line  about  the  optic  axis  must  be  adjusted.  The  simpler  system  in  Fig.  lb  is  preferred  for 
rapid  scanning  [15]  so  that  complex  alignments  do  not  complicate  the  delicate  high  speed  mirror 
drive(s).  With  this  design,  two  axes  of  tilt  for  one  mirror  must  be  aligned  with  respect  to  the 
other  so  that  wave  fronts  returning  from  both  mirrors  will  be  parallel  in  the  output  beam. 
However,  as  will  be  seen,  this  function  can  be  performed  by  auxiliary,  non-moving  mirrors.  Our 
design  for  a  rapid-scanning  interferometer  has  some  novel  features  but  is  closely  related  to  the 
design  shown  in  Fig.  lb. 

2.2  Our  Design 

The  arrangement  of  mirrors  and  beam  splitter  used  in  our  design  is  shown  in  Fig.  2.  The 
mirror  drive  is  somewhat  novel  for  this  type  of  interferometer  in  that  one  mirror  in  each  beam. 
M3  and  M4.  is  moved.  The  mirrors  move  in  opposite  directions,  since  they  are  mounted  on 
opposite  sides  of  a  parallelogram  formed  by  the  mirror  carriers  and  the  levers  LI  and  L2.  Each 
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lever  pivots  in  the  middle  with  respect  to  the  main  interferometer  drive  base,  which  is  shown. 
The  ends  of  the  levers  connect  to  the  mirror  carriers  also  through  pivots.  The  angle  of  L2  can  be 
changed  by  the  torque  motor  whose  armature  is  attached  to  the  lever  and  whose  field  is  attached 
to  the  base.  As  a  result,  the  scanning  mirrors  move  in  such  a  way  that  the  plane  of  each 
reflecting  surface  is  always  parallel  to  the  initial  plane,  and  the  path  difference  between  the 
beams  is  four  times  the  displacement  of  one  of  the  mirrors  from  the  position  where  they  are 
coplanar. 


Path  differences  up  to  8.2  cm  can  be  achieved  with  an  angular  rotation  of  about  20°.  This 
makes  a  fairly  compact  drive  for  such  a  resolution.  It  is  also  well-balanced  which  is  important  for 
keeping  vibration  down  in  rapid  scanning.  Mirror  speeds  up  to  5  cm/sec  can  be  obtained  with  the 
servo-control.  For  spectroscopy  near  sea  level,  absorption  lines  in  the  transmission  spectrum  of 
the  atmosphere  over  a  400  m  path  are  generally  of  a  such  a  width  that  7.5  GHz  resolution  is 
adequate.  This  corresponds  to  a  path  difference  of  4  cm  and  a  mirror  displacement  of  1  cm.  Thus, 
a  spectrum  can  be  obtained  in  around  200  msec. 


The  mechanism  of  this  drive  has  some  antecedents.  It  is  similar  to  a  design  for  translating  a 
single  mirror  in  an  interferometer  developed  at  AFGL  [9,10],  In  this,  only  one  mirror  was 
translated  by  levers  on  one  side  of  the  base.  The  mirror  also  had  a  shorter  range  of  travel  since 
the  pivots  were  a  tortional  type  limited  by  strain  failure.  We  tried  to  use  a  similar  pivot  (a  model 
2479167-Z  Flex  Pivot  made  by  Bendix,  Electric  Fluid  and  Power  Div.)  which  was  more  flexible 
and  able  to  reach  ±25°  with  indefinite  life,  provided  there  was  zero  load.  Such  pivots  were 
attractive  because  there  is  practically  no  friction  involved,  but  they  proved  too  flexible  in  two 
respects:  1)  they  could  not  hold  the  armature  in  the  center  of  the  stator  of  the  torque  motor,  and 
2)  the  scanning  motions  created  destructive  resonant  vibrations  in  the  flexible  pivots.  For  these 


and  other  reasons,  we  found  it  necessary  to  replace  the  pivots  with  jeweled  bearings  and  pivots, 
which  we  made  to  a  tolerance  of  about  8  pm  (0.0003  in).  A  light-weight  lubricant  (Synt-A-Lube, 
made  by  H.  Molbius  in  Switzerland)  was  used  for  acceptably  low  friction. 

As  may  be  seen  in  the  photographs  in  Fig.  3,  the  levers  LI  and  L2  connect  to  the  base  and 
to  the  mirror  carriers  with  a  top  and  bottom  bearing/pivot  set,  with  one  exception:  only  one 
bearing/pivot  is  used  to  connect  L2  and  the  mirror  carriers.  These  have  a  special  arrangement  to 
allow  for  an  important  adjustment  If  one  tries  to  assemble  the  scanning  mechanism  without 
making  allowance  for  possible  imperfections  in  machining  the  parts,  a  bearing/pivot  set  might 
not  fit  in  the  last  joint  or  might  be  seriously  strained.  To  alleviate  this,  it  is  necessary  to  have 
some  adjustment  in  the  length  of  at  least  one  of  the  levers  or  mirror  carriers.  However,  this  is  not 
necessarily  satisfactory  since  the  assembled  parts  may  not  then  form  a  parallelogram,  and 
parallel  motion  of  the  mirrors  may  be  lost.  Parallelograms  can  be  obtained  on  both  sides  of  the 
drive  if  an  adjustment  can  be  made  in  the  length  of  both  mirror  carriers  and  both  sides  of  the 
lever  L2.  The  adjustment  is  made  in  each  case  by  moving  the  hole  in  which  the  bearing  or  pivot 
fits  at  the  back  of  the  mirror  carriers  and  ends  of  L2.  This  hole  is  made  to  be  the  inner  hole  of 
two  nested  collets  which  are  each  slightly  and  equally  eccentric.  By  rotating  the  collets  with 
respect  to  the  mirror  carrier  or  lever,  one  can  move  the  center  line  of  the  hole.  When  a  position  is 
found  that  makes  the  mirrors  move  properly  and  which  relieves  strain,  the  collets  are  clamped. 
The  advantage  of  this  method  is  that  it  resists  loosening  during  scanning. 

Motion  of  the  mirrors  is  monitored  by  special  optical  grating  encoders  made  by  Teledyne 
Gurley.  The  gratings,  one  on  each  side,  are  mounted  to  the  interferometer  base  as  shown  in  Figs. 
2  and  3,  and  the  two  reading  heads  to  the  mirror  carriers.  The  gratings  have  an  envelope  which  is 
a  shallow  arc  corresponding  to  the  arc  of  the  levers.  (While  the  mirrors  move  so  that  the 
reflecting  plane  is  parallel  to  its  starting  position,  the  reader  will  notice  that  the  mirrors  move  in 
and  out  slightly  from  the  interferometer  base  due  to  the  arc  taken  by  the  ends  of  the  levers.  The 
mirrors  are  large  enough  that  they  always  reflect  the  radiation  beams  in  the  interferometer,  so 
this  is  not  a  problem.) 

The  period  of  the  grating  in  the  Teledyne  Gurley  encoders  is  20  pm,  and  with  the 
conditioning  electronics,  an  electrical  pulse  can  be  derived  for  every  5  pm  of  motion  of  the 
reading  head.  These  interpolated  pulses  are  triggered  by  the  zero-crossings  of  the  two  differential 
channels  in  quadrature.  The  accuracy  of  encoding  depends  essentially  on  two  factors:  (1)  the 
fundamental  accuracy  of  the  ruling  of  the  grating,  and  (2)  the  accuracy  with  which  mechanical 
and  electrical  adjustments  can  be  made  to  achieve  both  a  uniform  duty  cycle  between 
zero-crossings  and  the  phasing  between  the  channels  nominally  in  quadrature.  In  practice,  the 
latter  controls  the  accuracy  achievable,  and  the  manufacturer  estimates  that  this  is  about  ±2  pm. 
At  ±Xyi50  for  the  shortest  wavelength  of  possible  interest,  this  degree  of  accuracy  is  adequate. 
Encoder  gratings  are  attached  as  shown  in  Figs.  2  and  3,  and  adjustments  are  provided  for 
alignment  of  the  gratings  with  the  mirrors  and  with  the  the  reading  heads.  Measurements  from 
the  encoders  from  both  sides  of  the  drive  are  used  in  reckoning  path  difference. 

Folding  mirrors  Ml  and  M2  are  provided  so  that  the  two  moving  mirrors  can  be  parallel  to 
each  other  instead  of  perpendicular,  as  in  a  conventional  interferometer.  These  mirrors  also  make 
it  possible  to  align  the  interferometer  in  the  sense  of  bringing  the  output  beams  from  the  two 
moving  mirrors  into  coincidence.  This  is  tested  by  shining  a  helium  neon  laser  into  the  input  port 
and  checking  that  the  return  beams  are  in  coincidence  on  some  distant  target.  Diffraction  from 
the  grid  beam  splitter  is  not  a  problem  since  the  zero-order  spots  are  quite  distinct.  In  fact  the 
diffraction  patterns  are  useful  in  insuring  the  correct  orientation  of  the  grid  lines. 
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Fig.  3  Two  views  of  the  interferometer  scanning  mechanism.  Note:  two  bearings 
attach  the  levers  to  the  interferometer  drive  base  and  lever  LI  is 
connected  to  the  mirror  carrier  by  two  bearings.  Lever  L2  is  connected 
to  the  mirror  carrier  by  one  bearing  mounted  in  the  double  eccentric 
collets,  mentioned  in  the  text.  The  signal  conditioning  units  for  the 
encoders  have  been  removed  for  clarity. 


.*  A/./.  A.  A.  A.  A.  A.  A. A  A  A  A  ,  - 


3  OPTICAL  DESIGN 


3.1  Transmit/Receive  Layout  and  Polarization  Control 

An  input/output  polarizing  beam  splitter  is  an  important  ingredient,  in  addition  to  those 
shown  in  Fig.  2,  for  proper  operation  of  the  interferometer.  This  component,  shown  in  Figs.  4 
and  5,  has  three  functions:  (1)  it  establishes  an  input  plane  of  polarization,  (2)  it  acts  as  a 
polarization  analyzer  for  the  output  from  the  interferometer,  and  (3)  it  makes  it  possible  to  have  a 
monostatic  antenna  (for  transmitting  and  receiving  signals).  This  was  important  in  making  it 
easier  to  implement  a  short  path  for  calibration  (see  Section  6). 


The  signal  from  the  lamp  source  is  focused  and  reflected  up  to  the  input/output  beam 
splitter.  There,  one  plane  of  polarization,  at  45°  with  respect  to  a  vertical  line,  is  reflected  into  the 
interferometer.  The  beam  splitter  in  the  interferometer  has  vertical  grid  lines,  so  that  half  the 
incoming  radiation  is  split  to  each  of  the  sides  of  the  interferometer.  The  same  polarizations 
return  from  mirrors  M3  and  M4.  When  they  recombine  at  beam  splitter  BD2,  there  is  a  phase 
delay  between  the  beams  proportional  to  the  path  difference.  When  the  phase  difference  for  a 
given  wavelength  is  zero  or  a  multiple  of  271,  the  resultant  polarization  going  back  out  the 
interferometer  is  the  same  as  the  input.  For  phase  differences  which  are  odd  multiples  of  n,  the 
outgoing  polarization  is  perpendicular  to  that  in  the  input.  For  intermediate  phase  differences,  the 
polarization  is  elliptical  or  circular.  Regardless,  the  input/output  beam  splitter  acts  as  an 
analyzer:  it  passes  only  polarizations  perpendicular  to  that  from  the  input.  With  a  conventional 


interferometer,  the  output  signal  at  zero  path  difference  is  a  maximum.  For  the  interferometer 
described  here,  interferograms  are  inverted  but  this  has  only  a  trivial  effect  on  the  extraction  of 
spectra  by  Fourier  transformation. 

The  signal  leaving  the  input/output  beam  splitter  has  a  polarization  which  is  once  again  at 
an  angle  of  45°  with  respect  to  the  vertical.  This  signal,  which  passes  out  through  the  Cassegrain 
antenna,  propagates  from  one  building  to  another  on  the  Georgia  Tech  campus  over  a  path  200  m 
long.  At  the  end  of  this  path,  a  dihedral  retroreflector  has  been  sited.  Shown  in  Fig.  6,  its  roof 
edge  is  horizontal  so  that  the  plane  of  polarization  is  rotated  by  90°.  The  retroreflector  is  aimed 
back  at  the  transmitter/receiver  where  the  radiation  is  collected  by  the  Cassegrain  antenna.  The 
signal  returns  to  the  input/output  beam  splitter  where  it  is  now  reflected  upwards  to  the  detector. 
In  order  for  the  detector  to  see  this  signal  alone,  a  polarization  Filter  is  provided,  made  from  two 
polarizing  grids.  Otherwise,  the  detector  would  see  a  large  signal  directly  from  the  source.  While 
this  signal  is  unmodulated  by  the  interferometer,  it  might  overload  the  detector  if  not  removed. 

The  retroreflector  has  been  made  from  Hexcell  ™  honeycomb  aluminum  plates  which 
were  found  to  be  flat  to  a  few  thousandths  of  an  inch,  adequate  for  the  NMMW  region.  They 
were  mounted  in  a  machine  jig  to  insure  squareness  of  the  roof  intersection.  Proper 
retroreflection  has  been  verified  by  clamping  plate  glass  mirrors  to  the  retroreflector  surfaces. 
The  retroreflector  assembly  has  been  mounted  to  a  specially  constructed,  large  tripod  assembly, 
which  has  been  reinforced  with  guy  wires  and  shielded  from  wind  and  weather  by  a  scaffold  clad 
with  tarpaulin. 
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Fig.  6  Photograph  of  main  retroreflector 


3.2  Cassegrain  Optics 

The  Cassegrain  antenna  consists  of  a  630  mm  diameter  electroformed  primary  paraboloid, 
with  a  focal  ratio  f/0.28,  and  a  132  mm  diameter  secondary  hyperboloid.  A  photograph  is  given 
in  Fig.  7.  The  surface  of  the  mirror  is  nickel,  over-coated  with  rhodium,  chosen  for  weather 
resistance.  The  secondary  was  a  nickel-coated  aluminum  mirror  which  was  diamond  turned  by 
Aero-Research  Corporation.  This  combination  provided  optically  polished  surfaces  so  that  the 
antenna  could  be  aimed  at  the  remote  retroreflector  using  optical  images. 


Construction  of  this  type  of  antenna  requires  some  care  since  the  electroformed  primary  is 
relatively  thin  (0.08  in,  2  mm).  With  the  advice  of  the  maker  (Optical  Radiation  Corp.),  we 
mounted  the  primary  in  a  shell  that  supported  the  back  of  the  mirror  at  about  2/3  the  diameter  of 
the  mirror.  We  fastened  it  to  the  shell  using  a  specially  selected  silicone  adhesive  recommended 
as  likely  to  minimize  stress  on  the  mirror.  The  shell  in  turn  supported  the  secondary  mirror  with 
conventional  spars.  Since  it  was  impossible  to  guarantee  that  the  primary  would  be  located 
accurately  with  respect  to  the  secondary  holder,  the  latter  was  provided  with  adjustments,  both 
for  centering  and  for  tilt.  Subsequent  tests,  however,  showed  that  the  primary  had  been 
accurately  positioned  and  that  adjustments  from  the  neutral  position  of  the  secondary  were  not 
necessary. 

The  whole  antenna  subsystem  is  mounted  in  an  Aerotech  gimbal  with  micrometer 
adjustments  for  aiming.  We  built  a  pulley  counterbalance  system  to  stay  within  in  the  static 
torque  limitations  of  the  gimbal. 

3.3  Source 

The  source  of  radiation  used  in  the  system  is  a  hot  mercury  lamp  mounted  in  a  cavity  with 
a  conical  horn  (or  cone  channel)  output.  Such  lamps  have  been  found  to  radiate  in  the  NMMW 
band  like  a  black  body  of  2000K  [16].  The  cavity  is  127  mm  (5  in)  in  diameter  and  integrates  the 
signal  from  all  parts  of  the  lamp  into  an  apparent  source  about  5.3  mm  (0.25  in)  in  diameter  at 
the  output  port  of  the  cavity.  The  surface  of  the  brass  cavity  has  been  gold-plated  for  high 
reflectivity  at  NMM  wavelengths. 

The  copper  electroformed  cone  has  apertures  of  5.3  mm  and  43  mm  and  a  length  of 
1 16  mm.  It  serves  two  purposes:  1)  it  restricts  radiation  to  cone  angles  which  just  fill  the  optical 
train  after  the  source,  and  2)  it  makes  the  source  appear  to  be  larger  in  diameter  [17]  so  that  the 
image  of  the  source  projected  on  the  retroreflector  is  large  enough  to  fill  it.  We  have  estimated 
the  apparent  size  of  the  source  to  be  33  mm  in  diameter  and  the  image  at  the  remote 
retroreflector  to  be  3.5  m  in  diameter.  There  is  some  vignetting  expected  near  the  outer  part  of 
this  due  to  the  sizes  of  the  optical  elements,  but  the  central  1.6  m  of  the  image  should  be  free  of 
vignetting. 

As  shown  in  Figs.  4  and  5,  there  is  a  70  mm  diameter  relay  lens  after  the  cone  to  place  an 
intermediate  focus  of  the  source  on  the  input/output  beam  splitter  BD1.  The  lens  is  double 
convex  (hyperboloid)  and  made  of  TPX®,  a  plastic  selected  for  its  low  losses  in  the  NMMW 
band  and  for  its  transparency.  The  hyperboloid  surfaces  were  made  with  the  usual  eccentricity 
(equal  to  the  refractive  index)  [18]  to  minimize  spherical  aberration  for  on-axis  radiation. 

3.4  Interferometer  Optics 

The  collimating  lens  C2  in  the  interferometer  has  an  interesting  dual  role:  1)  it  forms  a 
nearly  parallel  beam  of  radiation  for  entry  into  the  interferometer,  and  2)  it  controls  the  quality  of 
the  quasi-collimation  by  the  Cassegrain  antenna.  It  is  a  plano-convex  (hyperboloid)  lens  made  of 
TPX®  and  is  100  mm  in  diameter. 

Collimation  to  the  interferometer  is  not  perfect  due  to  diffraction  and  to  the  extended  size 
of  the  source.  Since  the  optical  path  lengths  in  the  interferometer  are  not  large  and  since  the  ratio 
of  aperture  to  wavelength  is  generally  better  than  50:1,  diffraction  will  be  neglected  in  the 
interferometer.  The  mirrors  of  the  interferometer  are,  however,  somewhat  over-sized  to 
accommodate  a  small  amount  of  diffraction.  The  effect  of  the  extended  source  is  to  cause  the 


beam  envelope  in  the  interferometer  to  decrease  from  an  approximate  diameter  of  83  mm  at  C2 
to  a  diameter  of  45  mm  in  the  vicinity  of  the  moving  mirrors  M3  and  M4.  The  reason  for  this  is 
that  these  mirrors  were  placed  near  a  pupil  which  is  the  image  of  the  antenna  aperture  behind  C2. 

In  initial  designs  of  the  lens  C2,  the  surface  was  hyperboloid  with  the  usual  eccentricity  to 
minimize  axial  aberration  [18].  However,  it  is  known  that  off-axis,  such  aspheric  surfaces  can 
actually  cause  aberrations  to  be  worse  than  for  the  case  where  spherical  lens  surfaces  are  used 
[19].  In  our  optical  system  there  are  many  optical  elements  in  the  train  before  the  beam  which 
propagates  to  the  retroreflector.  To  estimate  the  effects  of  aberration,  we  used  computer  ray 
tracing.  We  attempted  to  focus  the  Cassegrain  so  that  there  was  an  image  formed  of  the  source 
on  the  retroreflector.  We  observed  that  the  image  field  was  drastically  curved:  the  position  of 
focus  for  the  apparent  edge  of  the  source  was  at  148  m  range,  as  compared  with  200  m,  on  axis. 
This  might  create  a  little  wash  out  of  the  image  near  the  edges  of  the  retroreflector  so  we  devised 
a  correction.  We  made  the  eccentricity  of  the  hyperboloid  somewhat  smaller,  until  by  trial  and 
error,  the  aberrations  were  negligible.  The  final  eccentricity  was  0.87  n,  instead  of  the  usual  n, 
where  n  is  the  refractive  index,  1.458,  for  TPX®.  It  is  interesting  that  it  was  possible  to  correct 
the  image  curvature  by  adjusting  the  surface  of  only  one  element  in  a  simple  way.  A  theoretical 
analysis  [20]  of  third-order  aberrations  in  the  system  showed  that  this  plano-convex  lens,  which 
radiation  passes  through  twice,  was  the  primary  cause  of  the  field  curvature,  and  the  special 
circumstances  of  the  layout  led  to  a  fortuitously  simple  solution. 

3.5  Detector  Optics 

The  optics  leading  to  the  detector  from  the  input/output  beam  splitter  BD1  consist  of  a 
36  mm  diameter  TPX®  lens  for  relaying  the  signal  from  the  retroreflector  and  two  polarizing 
filters  to  block  direct  radiation  from  the  source.  The  window  of  the  detector  is  13  mm  in 
diameter. 

4  DETECTION 
4.1  Detector  Type 

The  lowest  noise  detectors  sensitive  over  the  NMMW  range  at  the  time  this  work  began 
were  liquid  helium-cooled  composite  bolometers,  based  on  Germanium  [21].  Minimum 
detectable  power  values  (for  a  radiation  band  0-1000  GHz)  can  be  around  NEP=10  UW/Hz‘^. 
Regrettably,  these  detectors  cannot  be  used  in  this  application  because  their  response  rolls  off  at 
frequencies  above  (typically)  tens  of  Hz.  For  rapid  scanning  of  the  shortest  wavelengths  that 
might  be  seen  here,  the  interferometer  could  produce  modulations  of  signal  as  high  as  150  Hz  in 
frequency.  The  upper  limit  of  the  detector  modulation  frequency  response  should  be  significantly 
higher  to  avoid  phase  shifts.  Hot-electron  InSb  photoconductive  detectors,  on  the  other  hand, 
have  modulation  frequency  responses  up  to  at  least  500  kHz,  but  a  higher  noise  level.  For 
operation  at  liquid-helium  temperatures,  NEP  values  around  10U  W/Hzlrt  have  been  reported 
[22],  Values  of  this  order  can  be  shown  to  be  the  minimum  necessary  to  obtain  spectra  in 
reasonable  amounts  of  total  measurement  time,  as  discussed  in  Section  7. 

Obtaining  a  working  InSb  detector  with  such  sensitivity  proved  to  be  difficult  during  the 
contract  period.  Georgia  Tech  owned  InSb  elements  at  the  start  of  the  work  of  unknown 
sensitivity.  These  proved  to  be  significantly  less  sensitive  than  Golay  cells,  which  usually  have 
NEP  values  of  10’  to  lO10  W/Hz1/J.  We  then  were  grateful  to  obtain  an  element  through  MICOM 
and  learned  to  contact  it  to  the  post  which  holds  it  in  our  detector  cavity.  The  resulting  sensitivity 
was  better  but  not  adequate  to  see  the  source  signal  over  the  400  m  path.  After  recontacting  the 
element  several  times,  the  best  that  we  achieved  was  2.4x10  "  W/Hzlfl  (see  Section  8),  not 


sufficient  to  conduct  the  experiment.  We  attempted  to  improve  our  technique  by  contacting  the 
element  in  a  reducing  atmosphere  but  this  was  unsuccessful.  Nevertheless,  it  is  still  believed 
possible  to  obtain  a  detector  with  the  necessary  sensitivity,  particularly,  since  new  sources  for 
detector  elements  have  appeared  in  the  U.S.  since  the  work  began. 

4.2  Dewar  and  Internal  Optics 

The  detector  element  used  has  been  mounted  in  a  standard  Infrared  Laboratories  dewar 
with  a  downward  looking  window,  as  shown  in  Fig.  8.  The  detector  element  is  contacted  to  a 
post  within  a  gold-plated  beryllium-copper  cavity  which  is  fed  by  a  small  horn.  A  fine  wire  to  the 
other  side  of  the  element  completes  the  circuit  to  the  preamp. 


Dewar  cold  surface 


Outer  Dewar  jacket 


Fig.  8  Diagram  of  internal  detector  optics 


It  is  necessary  in  Fourier  spectroscopy  to  have  a  low-pass  optical  filter  to  prevent  aliasing 
of  high  frequency  signals  in  the  region  of  interest.  InSb  elements  have  a  natural  roll  off  at  around 
1000  GHz.  However,  it  is  necessary  to  remove  higher  frequencies  in  the  infrared  where  the 
element  can  be  sensitive.  Accordingly,  black  polyethylene  and  crystal  quartz  have  been  mounted 
on  the  radiation  shield  of  the  detector  dewar  and  on  the  detector  cavity.  The  external  detector 
window  is  made  of  PTFE. 

4.3  Preamp 

A  standard  low-noise  FET  preamplifier  with  three  stages  and  a  gain  of  300  has  been  used 
with  the  detector.  The  noise  figure  of  the  amplifier  was  tested  to  be  1.8  nv/Hzw.  This  compares 
favorably  with  3  nv/Hzlrt  typical  of  such  preamps. 
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5  ELECTRONIC  CONTROL 

The  electronics  subassemblies  prepared  for  the  FTS  are  diagrammed  in  Fig.  9.  Many  of  the 
electronic  systems  involved  were  developed  concurrently  for  this  program  and  for  research  for 
AFGL.  The  interested  reader  may  consult  [23,24]  for  detailed  information  on  the  circuits  and  on 
common  optical  elements.  Control  of  the  interferometer  is  separated  into  two  sections,  a 
primarily  analog  section,  called  "Head  Electronics",  near  the  interferometer,  and  a  primarily 
digital  section,  called  "Control  Electronics",  which  can  be  several  feet  away. 
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Fig.  9  Block  diagram  of  electronics  in  FTS 


Two  Motorola  6809  microprocessors  are  contained  in  the  Control  Electronics  chassis,  one 
responsible  for  servo-control  and  one  for  monitoring  keyboard  inputs  which  set  various  control 
parameters.  The  circuits  in  the  Head  Electronics,  in  conjunction  with  the  Control  Electronics, 
implement  a  hybrid  analog/digital  servo-system  for  accurate  control  of  the  moving  mirrors.  The 
primary  mode  of  operation  drives  the  mirrors  at  a  constant  speed  until  a  predetermined  path 
difference  is  reached  and  then  the  motion  is  reversed.  The  servo-control  circuitry  supplies  a 
control  voltage  to  a  power  servo-amplifier  used  as  a  voltage  controlled  current  source  (VCC) 
which  supplies  current  to  the  torque  motor.  The  mirror  position  is  determined  by  the  pair  of 
incremental  position  encoders,  whose  conditioning  electronics  closes  the  feedback  loop  by 
supplying  information  to  the  position  counter.  The  position  counter  supplies  position  information 
to  the  servo-loop  control  circuitry  and  to  the  control  processor. 

Fig.  10  gives  a  functional  block  diagram  of  the  servo-svstem  that  has  been  implemented.  In 
the  intial  design,  the  servo-loop  controller  was  directly  connected  to  the  servo-amplifier.  Late  in 
the  research  effort  it  was  discovered  that  the  digital  feedback  of  position  and  velocity 
implemented  by  this  system  caused  a  grainy  signal  to  be  transmitted  to  the  torque  motor.  As  a 
new  byte  of  position  or  velocity  information  was  received,  the  torque  motor  was  asked  to  make  a 
small  step  response.  This  let  the  motor  to  radiate  a  great  deal  of  interference  which  was  picked 
up  by  some  part  of  the  sensitive  detector  system.  As  a  result,  interferograms  had  an  interference 
component  superimposed,  as  shown  in  Fig.  11. 
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Fig.  10  Servo-system  functional  diagram.  In  the  upper  right  comer,  there  is  a 
switch  to  connect  either  the  servo-loop  controller  or  a  function 
generator  to  the  power  servo- amplifier.  In  the  initial  design  the 
servo-loop  controller  was  connected  and  later  designs  connected  the 
function  generator,  as  explained  in  the  text. 


Manual  rapid  scanning  showed  (see  Fig.  12)  that  this  spurious  signal  did  indeed  arise  from 
the  servo-system,  so  it  was  necessary  to  remove  the  continuous  digital  feedback.  A  centering 
spring  was  added  to  the  mirror  drive  and  a  function  generator  was  used  to  scan  the  interferometer 
in  an  open  loop.  This  mode  is  selected  in  Fig.  10  when  the  function  generator  is  switched  in. 

In  the  revised  system,  the  encoder  pulses  are  not  used  for  feedback  to  the  servo-system,  but 
they  are  used  by  the  control  processor  to  send  to  the  data  logging  system  to  trigger 
analog-to-digital  conversion.  This  is  the  primary  purpose  of  the  encoders  in  any  case.  With  the 
Control  Electronics,  one  can  define  the  size  of  the  sampling  interval  in  terms  of  the  number  of 
encoder  pulses  between  samples.  In  the  course  of  this  work,  various  data  logging  computers  have 
been  used  and  need  not  be  commented  on  in  detail.  The  same  computers  have  been  used  to 
implement  the  Fourier  transforms  [251. 


6  CALIBRATION 


To  determine  atmospheric  transmission  from  spectroscopic  measurements  of  transmitted 
signals,  it  is  necessary  to  know  the  signal  that  would  have  been  obtained  if  the  atmosphere  were 
not  present  and  absorbing.  The  source,  the  instrumentation,  the  geometry  of  the  propagation 
path,  and  atmospheric  attenuation  all  contribute  to  the  shape  and  magnitude  of  the  spectrum  of 
radiation  received  over  the  propagation  path.  In  a  laboratory  experiment  involving  an  absorption 
chamber,  it  is  often  possible  to  pump  out  the  atmosphere  and  record  a  spectrum  characteristic  of 
the  instrumental-type  effects.  The  atmospheric  transmission  is  then  the  ratio  of  spectra  with  and 
without  the  atmosphere.  Clearly,  we  cannot  do  this  in  the  open  atmosphere. 

In  monochromatic  experiments,  methods  have  been  devised  to  remove  most  instrumental 
effects  [26|.  Multiple  corner  cube  reflectors  are  set  up  at  different  ranges  from  the  transceiver. 
Ratios  of  signal  received  from  the  various  comer  reflectors  depend  on  the  differential 
atmospheric  attenuation  between  the  ranges  and  on  differential  diffraction  losses.  Since  the 
monochromatic  sources  can  be  approximated  by  point  sources,  the  differential  diffraction  losses 
are  calculable  and  the  signal  ratios  can  be  corrected  to  find  the  contribution  of  atmospheric 
transmission. 
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In  the  use  of  extended,  wide-band  sources,  diffraction  effects  are  not  practically  calculable. 
The  best  that  one  can  do  is  record  the  spectrum  of  signals  received  from  the  full  path  (400  m) 
and  another  over  a  short  path,  and  then  take  the  ratio.  The  resulting  spectral  ratio  is  an 
approximation  to  the  transmission  spectrum  with  two  things  which  are  unknown:  the  actual 
transmission  scale  at  a  given  frequency,  and  any  remaining  instrumental  trends  in  the  spectrum. 
The  root  source  of  these  uncertainties  are  questions  about  how  the  propagated  beam  diverges,  as 
a  function  of  frequency,  resulting  in  incomplete  capture  of  the  signal  by  the  retroreflector  and  the 
receiving  antenna.  On  the  other  hand,  such  ratios  of  long  and  short  path  spectra  should  remove 
any  narrow  instrumental  effects,  such  as  features  due  to  internal  interference  in  the  lamp 
envelope.  All  that  should  remain  of  frequency  dependent  instrumental  effects  are  slow  trends  in 
propagation  efficiency.  Since  an  extended  source  is  used  that  over-fills  both  retroreflector  and 
the  receiver,  even  these  should  be  minor  over  much  of  the  NMMW  band.  Thus,  it  should  be 
possible  to  obtain  spectra  which  will  show  whether  or  not  isolated  anomalous  atmospheric 
absorption  features  (the  third  characteristic  mentioned  in  the  introduction)  appear  between  the 
known  water  vapor  lines. 

By  making  measurements  during  different  atmospheric  conditions,  ratios  of  atmospheric 
transmission  can  also  be  obtained  with  an  absolute  scale.  Since  warm,  dry  weather  has  been  said 
to  be  less  likely  to  produce  anomalous  absorption,  such  conditions  can  be  used  as  a  reference  to 
test  hypotheses  about  the  existence  of  anomalous  phenomena. 

For  the  facility  described  here,  a  small  retroreflector,  like  that  at  the  remote  site,  has  been 
constructed  as  shown  in  Fig.  13.  It  can  be  rolled  into  place  and  then  keyed  to  the  antenna  for 
reproducible  location.  This  has  been  useful  also  for  testing  the  interferometer  and  detector. 
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7  ENERGY  BUDGET 


In  this  section,  we  will  calculate  what  is  expected  in  signal  levels  through  the 
instrumentation  and  over  the  propagation  path.  We  can  make  two  useful  estimates:  the  total 
signal  in  the  pass  band  of  the  system,  and  the  signal  in  one  resolution  width  (7.5  GHz).  Most  of 
the  signal  of  interest  is  between  60  and  450  GHz,  since  the  instrumentation  will  probably  have 
poor  response  and  throughput  below  this  range,  and  the  atmosphere  will  absorb  most  of  the 
signal  above  this  range.  To  facilitate  calculation  of  the  total  signal,  we  divide  the  whole  range 
into  three  intervals,  between  60-180,  180-310,  and  310-450  GHz.  The  boundaries  of  these 
intervals  are  natural  ones  determined  by  noticeable  water  lines.  The  signal  received  in  each 
interval  is  then  the  product  of 

•  integrated  power  from  the  source 

•  instrumental  transmission 

•  fraction  of  signal  received  by  the  retroreflector 

•  fraction  of  signal  received  by  the  receiving  antenna,  and 

•  mean  transmission  of  the  atmosphere  over  400  m. 

We  begin  with  the  signal  radiated  by  the  source.  The  Rayleigh-Jeans  approximation  for 
black  body  radiation  holds  in  this  range.  Thus,  the  power  AP  available  in  a  resolution  width  Av 
can  be  approximated  by 

APS  =  IkTAOsi1— 
c 


and  the  integrated  power  P§  between  two  frequencies  v,  and  v2,  in  Hz,  is  given  by 


=  2kTAQ 


(V2-V1) 

3c2 


where  k  is  Boltzmann’s  constant,  T  is  the  source  temperature,  A  is  the  area  of  the  source,  Q  is 
the  solid  angle  in  which  radiation  is  collected  from  the  source,  and  c  is  the  speed  of  light.  We 
have  used  the  following  values: 

7  =  2000# 


A  =  1.075*  10"  V 
G  =  0.064  lsrr 


Values  of  AP,  and  P  calculated  for  the  three  intervals  are  given  in  Table  1 . 


Table  1.  Factors  in  the  Expected  Signals  over  the  400  m  Propagation  Path 


Frequency  Interval,  GHz 

Factor 

60-180 

180-310 

310-450 

Source  power  in  a  resolution 
interval,  W 

4.6x10’ 

1.9x10’ 

4.6x10® 

Integrated  source  power  in 
each  interval,  W 

7.9x10-’ 

3.4x1 0-7 

8.7xl0-7 

Instrumental  transmission 

0.047 

0.047 

0.047 

Chopper  factor  when  total 
power  is  measured 

0.5 

0.5 

0.5 

Fraction  received  by 
retroreflector 

0.08 

0.08 

0.08 

Fraction  received  by  receiver 

0.045 

0.045 

0.045 

Mean  atmospheric 
transmission 

0.90 

0.77 

0.25 

SIGNAL  NEAR  THE 
CENTER  OF  EACH 

5.9x10 13 

2.5xl012 

1.7x10 12 

INTERVAL  (IN  A 
RESOLUTION  WIDTH, 

7.5  GHz),  W 

@  110  GHz 

@  245  GHz 

@  340  GHz 

TOTAL  SIGNAL 

RECEIVED  IN  EACH 
INTERVAL,  W 

6.0xlO12 

2.2x10“ 

1.8x10" 

TOTAL  ALL  INTERVALS, 
W 

4.6x10“ 

The  instrumental  transmission  represents  all  transmission  factors  within  the 
transmitter/receiver.  There  are  a  great  many  such  factors,  some  of  which  are  well  known  and 
some  which  are  uncertain.  Among  those  which  are  known  are 

•  a  factor  of  1/2  for  the  selection  of  one  plane  of  polarization 

•  a  factor  of  1/2  for  interferometric  modulation.  This  is  applied  even  when  total 
power  is  measured  since  this  type  of  interferometer  must  be  set  at  some  large  path 
difference:  in  this  case  only  about  half  the  frequencies  transmit  through  the 
interferometer.  When  the  interferometer  is  scanning,  a  given  wavelength  is 
transmitted  about  half  the  time. 

•  a  factor  of  1/2,  if  a  chopper  is  turned  on 

•  a  factor  of  approximately  0.75  for  reflective  losses  at  the  eight  lens  surfaces 
encountered  by  the  radiation. 

Among  those  which  are  not  well-known  are 
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•  the  cumulative  effect  of  imperfections  throughout  the  system  including  radiation 
lost  at  the  edges  of  apertures 

•  the  coupling  efficiency  to  the  detector  element 

•  the  interferometric  efficiency,  generally  found  to  be  in  the  range  0. 8-0.9. 

Thus  we  cannot  estimate  the  total  instrumental  transmission  with  any  degree  of  certainty.  Since 
the  unknown  factors  could  be  significant,  we  select  a  value  of  0.25  to  represent  all  unknown 
factors  here,  leaving  a  total  instrumental  transmission  of  0.047  (except  for  a  chopper  factor). 

The  fraction  of  the  power  received  by  the  retroreflector  can  be  estimated  as  the  ratio  of  the 
area  of  the  retroreflector  to  that  of  the  ideal  projected  image,  or  1  m2/(3.52  m)^0.08.  From  ray 
traces,  we  estimate  that  the  beam  leaving  the  retroreflector  will  expand  to  twice  the  size  of  the 
retroreflector  plus  the  size  of  the  transmitting  antenna  by  the  time  it  returns  to  the  receiver. 
Thus,  the  fraction  received  by  the  receiving  antenna  is  rc(0.315/2.63)2=0.045.  The  atmospheric 
transmission  values  are  rough  approximations  calculated  from  an  empirical  model  [27]  for  a 
relative  humidity  of  50%,  and  temperature,  288K. 

From  the  results  in  Table  1,  we  conclude  that  a  detector  sensitivity  (NEP)  of  about 
10I2W/Hz1/2  is  necessary  to  have  a  signal  to  noise  ratio  in  the  spectrum  of  about  100  in  an 
integration  time  of  about  1/2  hr,  or  a  signal  to  noise  ratio  of  50  in  an  integration  time  of  a  little 
over  1  sec.  The  best  detector  sensitivity  achieved  was  24  times  worse  than  these  figures.  With 
such  a  sensitivity,  the  time  to  reach  100:1  SNR  in  a  spectrum  would  have  been  unacceptably 
long,  nearly  300  hrs. 

8  RESULTS 

The  main  results  of  this  work  to  date  have  been  tests  of  the  equipment  developed.  A  few 
tests  have  already  been  mentioned  in  foregoing  sections;  e.g.,  tests  of  the  rapid  scanning 
servo-system  and  of  the  preamplifier.  In  this  section,  tests  of  the  antenna  system,  of  the  spectral 
response  of  the  FTS,  and  of  the  sensitivity  of  the  detector  will  be  described. 

After  the  antenna  was  assembled,  but  before  the  interferometer  was  fully  mounted,  we 
placed  a  projection  lamp  at  the  Cassegrain  focus  of  the  antenna.  In  nighttime  experiments,  we 
adjusted  the  antenna  and  focused  it  on  the  tarpaulin  of  the  remote  retroreflector  to  determine  the 
quality  of  the  focus.  While  all  the  optical  elements  are  polished  so  that  they  look  optically  shiny, 
they  are  not  figured  to  optical  wavelength  tolerances.  The  projection  lamp  had  a  filament  that 
was  about  3  mm  across,  so  if  the  antenna  were  optically  perfect,  the  spot  on  the  retroreflector 
would  have  been  about  0.32  m  across.  The  observed  spot  was  about  0.75  m  across,  suggesting 
that  the  size  of  the  blur  circle  of  the  antenna  is  about  0.45  m.  Such  excellent  performance  should 
add  little  to  the  losses  in  the  system. 

Once  the  interferometer  and  other  optics  were  mounted  on  the  table  attached  behind  the 
antenna,  it  was  no  longer  possible  to  use  this  method  to  aim  the  antenna  since  the  Cassegrain 
focus  is  near  the  beam  splitter  BD1.  Instead,  rough  initial  alignment  can  be  obtained  by  several 
steps:  First  we  look  for  the  crude  image  of  the  retroreflector  during  daylight.  If  this  proves 
difficult,  we  can  place  a  bright  lamp  in  the  middle  of  the  retroreflector  at  night  and  look  for  the 
image  of  this  source.  We  then  reverse  the  process  and  place  a  photo-flood  in  front  of  the 
Cassegrain  aimed  at  the  retroreflector  with  plate- glass  mirrors  covering  the  Hexcell™  plates.  The 
retroreflector  is  then  aimed  so  that  the  reflected  light  is  centered  on  the  hut  on  the  roof  housing 
the  transmitter/receiver.  Finally,  the  mercury  lamp  source  is  temporarily  replaced  with  a  90  GHz 
Gunn  source  which  sends  a  clear  signal  through  the  whole  system  and  permits  final  adjustment 
to  maximize  this  signal. 
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During  the  course  of  the  work,  interferograms  were  produced  to  verify  that  good  spectra 
could  be  obtained.  An  interferogram,  which  was  presented  earlier  in  Fig.  12,  was  obtained  from 
the  short  path  to  the  calibration  reflector.  The  corresponding  spectrum  is  shown  in  Fig.  14.  The 
locations  of  expected  water  vapor  lines  are  shown  and  account  for  most  of  the  features  in  the 
spectrum.  One  at  about  660  GHz  is  probably  an  example  of  features  due  to  interference  in  the 
detector  window.  Such  features  are  commonly  found  and  should  disappear  when  ratios  of  long 
and  short  path  spectra  are  taken. 


Throughoui  the  work,  we  battled  problems  with  detector  sensitivity.  Improvements 
gradually  came  through  better  contacting  and  making  a  better  thermal  contact  between  the 
detector  post  and  the  liquid-helium  temperature  detector  cavity.  During  the  latter  stages,  we 
devised  a  laboratory  test  for  detector  NEP.  This  involved  a  large  aperture,  temperature-controlled 
black  body  source  about  80K  higher  in  temperature  than  the  room.  We  measured  the  signal  from 
this  source  using  a  chopper  and  a  lock-in  amplifier.  We  estimated  the  limiting  area-solid  angle 
product  in  the  system  and  calculated  the  optical  NEP  to  be  2.4X10-'1  W/Hz1/J.  The  voltage  noise 
of  the  cold  detector  element  was  measured  to  be  2.1  nv/Hzw.  Thus  the  responsivity  of  the 
detector  is  1 17  v/W. 

The  IV  characteristics  of  the  detector  were  also  measured.  The  dynamic  impedance  of  the 
detector  at  the  operating  current  of  70  microamps  was  found  to  be  272  ohms.  The  electrical  NEP 
derived  from  the  shape  of  the  IV  curve  was  found  to  be  5xl(>13  W/Hzlft.  This  value  is  in  line  with 
typical  InSb  detector  performance. 

The  size  of  the  discrepancy  between  the  optical  and  electrical  NEP  values  is  indicative  of 
three  possible  problems: 

•  pathological  low  sensitivity  to  radiation  due  to  poor  absorption  by  the  InSb 
element, 

•  a  bad  electrical  contact 

•  the  element  being  starved  for  radiation  in  its  particular  cavity  environment 

The  last  in  this  list  was  considered  the  least  likely,  and  it  is  difficult  to  test  for  the  first,  so  we 
attempted  to  improve  the  contacting.  It  is  known  that  the  best  InSb  contacts  are  made  in  an 
atmosphere  free  of  oxygen  and  water  vapor.  Thus,  we  attempted  to  set  up  the  entire  contacting 
process  in  a  glove  box.  However,  repeated  processing  degraded  the  element  to  the  point  where  it 
was  no  longer  useable.  Time  and  financial  resources  ran  out  in  this  contract  effort  before  this 
problem  could  be  resolved.  At  this  point,  the  best  course  of  action  will  be  to  purchase  a  new 
commercial  element  with  a  matching  cavity  and  feed. 

9  CONCLUSIONS 

In  the  course  of  this  work,  a  number  of  improvements  have  been  made  in  instrumentation 
for  NMMW  research.  In  particular,  an  interferometer  has  been  made  with  the  following 
combination  of  features: 

•  High  interferometric  efficiency,  due  to  the  Martin-Puplett  type  design 

•  Ruggedness  in  photolithographically  produced  beam  splitters 

•  Compactness  and  light-weight,  of  importance  in  attaching  it  to  an  antenna 

•  Capability  for  rapid  scanning,  enabling  the  recording  of  spectTa  in  a  fraction  of  a 
second. 

Admittedly,  severe  problems  have  occurred  in  this  work,  of  which  the  pivotal  one  in  the  end  was 
detector  sensitivity,  as  detailed  above.  Also,  much  more  time  had  to  be  spent  on  developing  the 
microprocessor-based  controls  for  the  system  than  had  been  anticipated.  As  a  result,  the  intended 
atmospheric  transmission  research  could  not  be  carried  out  inside  the  resources  of  the  contract. 
The  research  questions  that  were  raised  in  the  proposed  research  have  been  partially  answered 
since  the  inception  of  this  work.  Subsequent,  atmospheric  transmission  or  emission 
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measurements  have  largely  refuted  the  existence  of  anomalous  absorption  [e.g.,  6,7,27,28].  It 

would  still  be  desirable  to  see  wide-band  spectra  for  the  perspective  this  would  give,  and  it  is  our 

belief  that  the  any  remaining  problems  with  the  FTS  equipment  described  here  can  be  overcome. 
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Appendix  A 


Fabrication  of  Polarizing  Beam  Splitters 

by 

A.  Me  Sweeney 

The  polarizing  beam  splitters  consist  of  a  pattern  of  parallel 
strips  and  spaces  etched  in  an  aluminum  film  on  a  Mylar*  substrate.  The 
widths  of  the  opaque  strips  and  of  the  spaces  are  about  5  i^m  each. 
Therefore  there  are  1000  line-pairs  per  cm.  The  thickness  of  the 
aluminum  film  is  nominally  0.4  pm.  The  thickness  of  the  Type  S  Mylar* 
is  6  urn. 

The  general  procedure  for  fabricating  the  beam  splitters  is  as 
fol 1 ows : 

1.  A  piece  of  Mylar  is  stretched  over  an  aluminum  beam  splitter 

ring. 

2.  An  aluminum  film  is  deposited  on  the  Mylar  by  evaporating 
aluminum  wire  in  a  vacuum  chamber. 

3.  A  thin  layer  of  photoresist  is  applied  over  the  aluminum  layer. 

4.  The  photoresist  is  exposed  to  UV  radiation  through  a  mask. 

5.  Developing  the  resist  leaves  the  line  pattern  of  the  mask  in 
the  resist.  The  pattern  consists  of  parallel  strips  of  resist  separated 
by  uncoated  strips  of  aluminum. 

6.  The  strips  of  uncoated  aluminum  are  etched  away  from  the  Mylar 
leaving  the  pattern  of  opaque  aluminum  strips  separated  by  transparent 
spaces. 

♦Mylar  is  a  registered  trademark  of  DuPont. 


The  detailed  procedure  for  fabricating  the  beam  splitters  is  as 
f ol  1 ows : 

1.  A  piece  of  tylar  Is  stretched  over  an  aluminum  beam  splitter 
ring.  Installation  of  a  clamping  ring  further  stretches  the  Mylar. 

2.  A  measured  quantity  of  high  purity  aluminum  wire  is  placed  on  a 
tungsten  filament  in  a  vacuum  chamber.  The  stretched  Mylar  film  is 
cleaned  with  methanol  and  baked  dry.  It  is  then  placed  in  the  vacuum 
chamber  about  15  cm  above  the  filament  loaded  with  aluminum  wire. 

After  the  chamber  is  evacuated,  current  is  applied  to  the  filament 
to  first  melt  the  aluminum  wire  and  then  to  evaporate  it.  Some  of  the 
aluminum  vapor  coats  the  ffylar.  The  filament  is  kept  at  a  high 
temperature  until  most  of  the  aluminum  has  evaporated  from  the  filament. 

3.  A  resist  promoter,  hexamethyl disilazane  (HMDS)  is  then  applied 
and  spun  off  in  a  temperature  and  humidity  controlled  chamber. 

4.  Thinned  Shipley  AZ1350J  Positive  Resist  is  then  deposited  on 
the  aluminum  film.  The  resist  mixture  consists  of  two  parts  resist  to 
one  part  AZ  Thinner.  When  the  aluminum  surface  is  covered  with  the 
resist  solution  the  beam  splitter  ring  is  spun  at  about  3000  rpm  for  25 
seconds.  This  is  done  in  a  dust- free  environment.  The  resist  is  then 
"soft-baked"  at  85°C  for  25  minutes. 

5.  The  mask  and  resist  surface  are  then  pressed  together  in  an 
evacuable  container  shown  in  Figure  A1 .  The  container  has  a  large 
aperture  In  the  metal  bottom  to  allow  exposure  of  the  photoresist  to  UV 
radiation.  A  12  mm  inch  thick  transparent  fused  quartz  plate  is  placed 
over  an  0-ring  to  serve  as  a  window.  The  cleaned  mask  is  then  placed  on 
top  of  the  fused  quartz  plate  with  the  chromium  mask  pattern  facing  up. 


Then  the  beam  splitter  ring,  with  the  resist  surface  facing  down,  is 
placed  on  top  of  the  mask.  A  lapped  aluminum  disk  is  then  placed  inside 
the  beam  splitter  ring,  with  the  lapped  surface  pressing  down  on  the 
Mylar.  The  thickness  of  the  aluminum  disk  is  greater  than  the  thickness 
of  the  beam  splitter  ring,  so  the  top  of  the  disk  is  the  tallest  part  of 
the  stack.  The  outer  diameter  of  the  disk  is  just  slightly  smaller  than 
the  inner  diameter  of  the  beam  splitter  ring.  Therefore,  the  entire 
usable  area  of  the  photoresist  surface  is  pressed  into  contact  with  the 
mask.  Finally,  a  sheet  of  rubber  is  clamped  by  the  edge  to  the  top  of 
the  evaculable  container.  When  the  air  is  pumped  out  of  the  container, 
the  rubber  sheet  is  pressed  in  by  atmospheric  pressure.  It  presses 
against  the  aluminum  disk  which  in  turn  presses  the  flexible  Mylar  with 
its  aluminum  and  photoresist  coatings  into  contact  with  the  mask.  A  UV 
lamp  is  then  used  to  expose  the  photoresist  through  the  window  in  the 
bottom  of  the  container  and  through  the  mask  pattern.  The  exposure 
lasts  about  1  minute. 

6.  The  beam  splitter  ring  is  then  removed  from  the  evacuable 
container  and  the  photoresist  is  developed  in  a  mixture  of  20  parts  AZ- 
351  Developer  and  70  parts  water  for  15-20  seconds.  After  a  deionized 
water  rinse,  it  is  dried  and  examined  under  a  microscope. 

7.  The  photoresist  surface  is  then  "hard-baked"  at  110°C  for  20-30 
minutes. 

8.  The  aluminum  etch  solution  is  a  mixture  of  one  part  Acetic 
Acid,  16  parts  Phosphoric  Acid,  one  part  Nitric  Acid,  and  two  parts 
water.  A  drop  of  Kodak  Foto  Flo  is  added  as  a  wetting  agent  and  the 
temperature  is  raised  to  42.5°C.  It  takes  about  2  minutes  for  the 
aluminum  to  be  etched  away  sufficiently  to  see  light  through  the  beam 
spl  itter. 


9.  The  remaining  photoresist  is  removed  from  the  aluminum.  This 
is  accomplished  by  the  use  of  acetone.  The  beamsplitter  is  then  rinsed 
in  deionized  water  and  blown  dry  with  nitrogen. 


10.  The  final  step  is  to  remove  the  Mylar  film  from  the  aluminum 
ring,  clean  the  aluminum  and  the  Q-ring,  and  then  re- stretch  the  Mylar 
in  the  aluminum  ring. 


Appendix  B. 


THEORETICAL  OPTICAL  PROPERTIES  OF  A  GRID  ON  A  DIELECTRIC  FILM 

R.  A.  Bohlander 

1.0  Introduction 

A  grid  of  metal  strips,  or  wires,  preferentially  reflects  the  plane 
of  polarization  parallel  to  the  strips  and  preferentially  transmits  the 
orthogonal  polarization  when  the  wavelength  is  long  compared  with  the 
grid  period.  Polarizers  and  polarization  sensitive  beam  splitters  have 
been  constructed  on  this  principle  mainly  for  use  at  infrared  and  longer 
wavelengths.  When  the  application  is  in  the  microwave  region,  strips  or 
wires  can  be  stretched  on  a  frame  to  make  a  "free  standing"  grid.  This 
technique  has  been  extended  well  into  the  submillimetre  or  far  infrared 
region  (Ade  et  al.  1979;  Simonis  1979)  by  the  fabrication  of  grids  with 
periods  as  small  as  5  pm  and  wire  diametersas  smallas  12.5  pm. 
Transmissions  and  reflections  of  the  respective  polarizations  in  excess 
of  95  percent  have  been  measured  (Ade  et  al.  1979)  for  wavenumbers  as 
high  as  100  cm  However,  there  are  potential  advantages  in  forming 

the  grid  on  a  dielectric  substrate.  This  can  be  done,  for  example,  by 
photoetching  a  grid  pattern  in  a  thin  metal  layer  deposited  on  the 
dielectric.  The  substrate  then  helps  to  preserve  the  accuracy  of  the 
grid,  particularly  the  fine  ones  used  at  high  wavenumbers. 

The  performance  of  such  grids  has  been  studied  by  Auton  (1967)  in 
the  wavelength  range  20  to  100  pm  and  for  grid  periods  of  -1  to  10  urn. 
Each  disc,  made  of  polyethylene,  was  thick  relative  to  the  wavelength, 
and  its  surfaces  were  not  sufficiently  parallel  to  cause  resonances 
between  multiply  reflected  waves.  In  this  case,  the  dielectric  has  a 
deleterious  effect  on  transmission  as  summarized  by  Auton  (1967): 
"supporting  a  grid  on  a  substrate  of  refractive  index  n  not  only  reduces 
the  wanted  polarization,  as  compared  with  the  unsupported  grid,  but  also 
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increases  the  transmission  of  the  unwanted  polarization".  However,  if 
the  faces  of  the  substrate  are  parallel,  interference  between  multiply 
reflected  beams  can  occur.  If  the  thickness  of  the  substrate  is  chosen 
so  that  constructive  interference  gives  a  maximum  transmission,  the 
effect  of  the  substrate  on  the  optical  properties  of  the  polarizer  can 
be  negligible.  Auton  (1967)  commented  on  this  but  did  not  give  a 
theoretical  treatment. 

The  benefit  of  constructive  interference  can  also  be  realized  over 
a  wide  spectral  region  in  which  the  wavelength  is  much  longer  than  the 
thickness  of  the  substrate.  One  region  of  interest  is  the  near¬ 
millimetre  to  submillimetre  wavelength  range  where  new  component 
development  is  being  actively  pursued.  There  is  the  need  for  polarizing 
beam  splitters,  for  example  in  polarization  diplexers  (Chu  et  al.  1975) 
in  Fabry-Perot  etalons  (Ulrich  et  al.  1963),  and  in  the  Ma rtin-Puplett 
variation  (Martin  and  Puplett  1969,  Challener  et  al.  '.980)  of  the 
Michelson  interferometer.  Films  which  are  sufficiently  thin,  and  have 
sufficient  mechanical  strength  and  transparency  are  available 
commercially;  for  example:  Mylar  (Smith  and  Lowenstein  1975),  thickness 
i  2  urn?  Kapton  (Smith  and  Lowenstein  1975),  ^8  ym;  and  polypropylene 
(Chantry  and  Chamberlain  1972),  ^12  ym.  It  may  be  mentioned  that  these 
films  have  also  been  important  as  si±>strates  for  metal  meshes  used  in 
the  fabrication  of  multilayer  interference  filters  for  the  far  infrared 
(Holah  et  al.  1979). 

There  has  been  considerable  interest  in  the  theory  of  the  optical 
properties  of  grids.  Larsen  (1962)  reviews  work  up  to  about  1962  mainly 
on  grids  of  wires.  Some  of  the  more  recent  work  is  referenced  by  Kalhor 
(1978).  When  one  is  undertaking  the  design  or  the  evaluation  of  a 
polarizing  beam  splitter,  a  complete  picture  of  both  its  transmission 
and  reflection  properties  is  needed,  and  in  many  cases  arbitrary  input 
polarization  and  incidence  must  be  considered.  We  have  found  that  a 
more  complete  theoretical  model  is  needed  for  a  grid  of  metal  strips  on 
a  dielectric  film  than  is  readily  available  in  the  literature. 
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The  preliminary  account  which  is  given  here  is  an  extension  of  a 
theoretical  treatment  given  by  Adonina  and  Shestopalov  (1963),  hereafter 
referred  to  as  AS.  The  present  results  are  part  of  a  larger  effort, 
including  measurements,  which  will  be  reported  elsewhere  in  due  course. 

2.0  Theory 

2.1  No  Dielectric  Substrate 

We  will  consider  a  grid  which  has  a  period  g  and  metal  strips  of 

width  a,  as  shown  in  Figure  1,  and  initially  we  assume  that  the  metal 

strips  are  very  thin,  ideal  conductors,  Agronovitch  et  al.  (1962)  were 
able  to  express  the  diffraction  from  a  free  standing  grid  of  metal 
strips  as  a  boundary-value  problem  of  the  Riemann-Gilbert  type  for  which 
solutions  ware  known.  Their  approach  can  in  principle  give  a  full 
description  of  all  diffracted  orders.  The  validity  of  their  method  has 
been  checked  by  Kalhor  (1978)  and  compared  with  an  exact  solution  for  a 

special  case  by  Baldwin  and  Heins  (1954).  The  work  of  Agronovitch  et 

al.  was  a  seminal  paper  in  the  Soviet  literature  on  diffraction  by 
optical  systems  involving  metal  strips,  and  the  paper  by  AS  was  an 
extension  of  their  theory  to  the  case  of  a  dielectric  substrate. 

However,  polarizing  beamsplitters  are  used  at  wavelengths  which  are 
long  enough  that  only  zero  order  diffraction  need  be  considered.  We 
will  assume  that  the  wavelength  X  of  the  radiation  is  much  greater  than 
g  and  will  discuss  later  the  theoretical  limit  to  the  validity  of  this 
assumption.  It  is  helpful  to  subdivide  the  problem  into  cases  where 

a)  the  plane  of  incidence  is  either  perpendicular 
or  parallel  to  the  strips 

b)  the  plane  of  polarization  is  either  perpendicular 
or  parallel  to  the  strips. 

An  arbitrary  wave  can  be  subdivided  into  waves  of  the  two  polarizations 
mentioned  in  (b),  but  a  subdivision  for  an  abritrary  direction  of  incidence 
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along  the  lines,  in  (a)  is  of  course  not  possible.  Nevertheless,  it  is 
helpful  to  focus  on  the  two  limiting  cases  mentioned  in  (a)  for  the 
present.  For  the  amplitude  reflection  and  transmission  coefficients  r 
and  t  to  be  discussed,  we  will  use  a  first  subscript  to  specify  the 
orientation  of  the  plane  of  incidence  ,  and  the  second  ,  the  orientation 
of  the  polarization;  a  1  is  used  for  perpendicular  and  2  for  parallel. 

When  there  is  no  dielectric  substrate,  the  amplitude  reflection  and 
transmission  coefficients  in  the  long  wavelength  limit  are  well-known. 
Lamb  (1898),  Auton  (1967)  and  others  give  expressions  for  normal 
incidence  and  Marcuvitz  (1951),  Smirnov  (1958)  and  others  give  relevant 
expressions  for  oblique  incidence.  This  can  be  summarized  by  the 
following 


i  7  *  - 

l+i£  lnfu+l\cos6 
X  2 

Sl 

t  «  -i  A  cosQ 
12  l+ij>  ln(ofl)cos9 
X  2 


where  u=cos  [  tt  (  g  -  a)/g]and0is  the  angle  of  incidence.  Since  g  <<:*., 
reflection  dominates  for  this  polarization,  as  expected.  The 

coefficients  for  the  other  plane  of  incidence  r^^  and  t^  are 
obtained  from  and  t^  by  setting  the  factors  cos  6  to  unity.  The  case 

of  opposite  polarization,  ie.  r^,  t^  ,  etc.  may  be  obtained  from  the 
first  by  Babinet's  principle  (Born  and  Wolf  1970,  pp.  559-60;  AS).  For 
example,  r^  is  obtained  from  t^  but  with  u  replaced  by 
v  =  cos  (t:  a/g).  Intensity  reflection  and  transmission  coefficients  are 
given  by  rr*  and  tt*,  where  *  denotes  the  complex  conjugate,  and  these 
are  the  quantities  usually  tested  by  experiment.  Excellent  agreement 
between  experiment  and  theory  has  been  found  (AS;  Pursley  1956). 
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2.2  Dielectric  Substrate.  Only 


Again,  the  theory  is  well-known  (Born  and  Wolf  1970,  pp.  60-62). 
The  reflection  and  transmission  coefficients  for  a  dielectric  slab  which 
has  parallel  faces,  and  a  thickness  h  are 


r  ,r  (2iy-h) 

r  =  rg+  oexp  2  (2a) 

l+r^exp(2iY2^) 

t  t'  (iY0h) 

t  ,  a  oexp  2 _ 

S  l+r^exp(2iY,h)  ^b) 

a  t 

where  r„  is  the  amplitude  reflection  coefficient  of  each  surface  and  tj 
and  t'  are  the  amplitude  transmission  coefficients  of  the  first  and 
second  surfaces.  When  the  plane  of  polarization  is  perpendicular  to  the 
plane  of  incidence, 

r  -  y1~y2 

°  v4  (3a) 


w 


Yl+Y2' 


where  «  k  cos  9,  y2  *  kn  cos  6',  k  »  2  *  /  *  ,  n  is  the  index  of 

refraction  of  the  slab  (assumed  to  be  in  air),  and  9'  is  the  angle  of 
refraction.  For  the  opposite  polarization,  n  cos9'  is  replaced  by  cos 
in  equations  (3).  If  the  medium  is  absorbing,  n  is  replaced  by  the 
complex  refractive  index  n  +  iK  a/(4tt  )(  where  K  is  the  absorption 
coefficient.  Cos  6'  is  also  complex  and  n  cos  9'  may  be  expressed  as 
(Born  and  Wolf  1970,  pp.  627-30) 
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where  „  2  2  2  2  /  ,  2  2  2  2  2  2 

2v  =  n  -  »'  -  sin  6  +  /  (n  -<  -sin  6)  +  4>:  n 


2  ,2  2  ■  2qN  ^  /,  2  2  .  2  N 2  x  ,  2  2 

;  =  - (n  —  k  -  sin  9)  +  /  (n  -r  -sin  6)  +  4<  n 


k  *  KX/ (4ir)  . 


The  following  familiar  optical  properties  of  a  dielectric 
slab  carry  over  to  some  extent  when  a  grid  is  deposited  on  it: 

a)  as  the  thickness  increases,  the  reflection  and  transmission 
coefficients  cycle  between  maxima  and  minima  due  to  the 
interference  of  multiply  reflected  waves, 

b)  reflection  and  transmission  always  depend  on  the  angle  of 
incidence  (which  was  not  always  the  case  for  a  free 
standing  grid) 


at  Brewster's  angle  of  incidence,  6  ■  tan  ( 
the  dielectric  reflection  coefficient  is  small 
when  the  plane  of  polarization  is  parallel  to 
the  plane  of  incidence. 


The  last  property  can  be  used  to  enhance  the  degree  of  polarization 
obtained  from  a  grid  formed  on  a  substrate. 

2.3  Metal-strip  Grid  on  a  Dielectric  Slab 

It  is  desirable  to  have  expressions  for  the  optical  properties  of 
such  a  grid  for  arbitrary  planes  of  incidence.  Wait  (1957a,  b,  1959; 
see  also  Wait  1954,  and  Larsen  1962b)  has  done  this  for  grids  made  from 
wires.  For  metal  strips  on  a  dielectric  slab,  AS  have  considered  only 
the  plane  of  polarization  perpendicular  to  the  strips.  This  orientation 
gives  the  best  polarization  properties  at  oblique  incidence  for  two 
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reasons:  (1)  as  0  increases,  the  grid  period  appears  to  get  smaller 

relative  to  A  [(notice  the  cos  0  factors  in  eq.  (1)]  and  the  leakage 

of  the  "wrong  polarization"  decreases;  (2)  the  transmission  of  the 
"right  polarization",  namely  that  perpendicular  to  the  strips,  can  be 
improved  if  0  is  selected  to  be  Brewster's  angle  for  which  the 
dielectric  reflection  is  negligible.  The  latter  was  also  pointed  out 
recently  by  Challener  et  al.  (1980).  We  will  give  a  more  complete 
formulation  than  did  AS,  but  for  the  present  we,  too,  will  restrict  our 
attention  to  the  case  where  the  plane  of  incidence  is  perpendicular  to 
the  strips.  The  formulation  is  also  restricted  in  another  way;  namely, 
the  incident  beam  is  assumed  to  strike  first  the  side  with  the  metal  strips. 
Of  course,  the  intensity  transmission  coefficients  in  both  directions 
must  be  the  same,  on  thermodynamic  grounds,  but  the  reflection  and 
absorption  properties  need  not  be.  Other  directions  of  incidence  will 
be  covered  in  a  subsequent  paper. 

2.3.1  Polarization  Parallel  to  the  Strips 

Formulas  derived  by  AS  for  the  polarization  parallel  to  the  strips 
are  the  following: 

1  £  ln(u+l)  (Y-,-X0Y7)  -  k 

r12  =_2_J _ 2  _  ,  (5a)* 

i.  £.  ln(u+l)  (y1+X0Y2^  +  k 

2  A  2 


where 


=  Y2+Y1  -  (y2”Yj)  exp  (2iY2h) 
Y24*rl  +  ^Y2~Y1^  exp  (2i>2h) 


r12  +  1 


where  G  depends  only  on  the  dielectric  properties  and  is  given  by 


G-i 

0 


V^l  ^i(v:-v1)h^Y2~Yl  ei(Y2+Yl)h 
Y?  6  Yo  6 


*Note:  the  reader  should  be  warned  that  there  are  many  errors  in  the 

formulas  of  AS  which  appear  to  have  occurred  in  type-setting.  Corrected 
formulas  are  indicated  with  an  asterisk. 


When  there  is  no  dielectric  (ie.,  in  the  limit  when  n  1), 
equation  (5)  can  be  reduced  to  equation  (1).  Alternatively  when  there 
is  no  grid  (ie.,  in  the  limit  a  -*■  o) ,  equation  (5)  can  be  reduced  to 
equation  (2).  Also,  the  intensity  coefficients  r^2rl2*  anc *  ti2t12* 
can  be  shown  to  be  complementary;  ie., 


r  r  * 
12  12 


1  —  t  t  * 
12  12  ‘ 


An  assessment  of  the  limit  to  the  validity  of  the  long  wavelength 

approximation  can  be  made  by  comparison  of  equation  (5)  with  a  more 

complete  calculation  involving  higher  orders  of  diffraction  made  by  AS 

for  the  case  where  there  is  no  dielectric  slab.  The  comparison  is  shown 

in  Figure  2  for  two  angles  of  incidence.  When  the  incidence  is  near 

normal,  the  minimum  wavelength  X  where  there  is  good  agreement  is 

min 

X  .  -  2.5g.  As  9  increases,  X  .  increases  so  that,  for  9  =  65  , 

min  °  min 

X  min =  5g .  Unfortunately,  AS  have  not  provided  calculated  transmission 

values  in  the  case  where  the  grid  is  on  a  dielectric  slab.  However,  it 

is  reasonable  to  suppose  that  X  .  specified  in  air  should  increase 

min 

with  increasing  n  since  the  wavelength  in  the  dielectric  medium  is 
reduced . 


2.3.2  Polarization  Perpendicular  to  the  Strips 

The  formula  given  by  AS  for  reflection  is  as  follows 


if  ln(l+v)Vk n 

if  n2+l  lnO+v!Vkn+k 
X  2  2 


1,1  2  v,  R  +1 

—  +  n  1  o 
"  2 


Vn  Y1  exp(2iY,h) 

V.2t1 


‘V  .7 


AS  do  not  give  a  formula  for  t^,  but  by  analogy  with  equation  (5c)  it 
should  be  of  the  form 


(6d) 


In  the  limit  n  -*■  1,  G'  -*■  1  and  equation  (6d)  conforms  to  Babinet's 
principle;  in  the  limit  that  the  grid  is  absent  (a  -*■  0),  equation  (6d) 
must  go  over  to  equation  (2b).  This  fact  provides  a  means  of  deriving  a 
formula  for  G'  which  is 


2  2  2  2 

1  (n  Yj+y2)  ~  (n  Y1  -  y2)  exp  (2iY?h)  (6e) 


4  >i  y  2  n  exp  (i  y  2b) 


As  an  example,  we  show  in  Figure  3  calculations  of  the 
reflection,  transmission  and  absorption  of  a  grid  with  the  following 
pa rameters 

g  *  20  um 

a  =  10  um 

h  =  13  um 

n  *  1.5  (polypropylene) 

K  *  given  by  Chantry  and  Chamberlain  (1972) 

9  =  30°. 

2.4  Absorption 

Absorption  in  the  dielectric  substrate  may  be  included  in  equations 
(5)  and  (6)  by  following  the  procedure  for  replacing  n  and  n  cos  #  by 
their  complex  values  as  described  in  section  2.2.  However,  absorption 
by  the  metal  strips  has  not  as  yet  been  included  in  the  theory. 


properties  of  polarizer  (see  caption) 


Wavenumber  (cm  ) 

Figure  B. 3  Theoretical  optical  properties  of  a  metal  grid  polarizer  on  a 

polypropylene  film.  [Parameters  assumed:  g  =  20  urn,  a  =  10  _n 
h  *  13  pm,  0  =  30°,  n  *  1.5,  K  in  the  range  0  to  0.36cm  ^ 
(Chantry  and  Chamberlain  1972),  and  a  *  1  x  10^7  statmho/cn. 
Curve  a.)  rnru*  and  curve  b.)  1-r^/  -  tutu* 

curve  c.)  ohmic  losses  A. curve  d.)  t,  .t  *  and  1-r, „r  *• 


curve  e.)  l-r12r12  -  t12t12  ] 


12  12 


12  12 


M.M.' 


Generally,  this  is  a  small  effect  unless  the  beam  splitter  is  part  of  a 
strongly  resonant  cavity  such  as  a  Fabry-Perot  etalon.  According  to 
Ulrich  et  al.  (1963),  the  fraction  of  the  intensity  absorbed  by  the 
strips  is  of  the  order  of 


for  the  polarization  indicated,  where  o  is  the  conductivity  of  the 
strips  and  c  the  speed  of  light.  For  the  opposite  polarization,  ohmic 
losses  are  considered  negligible  (Ulrich  et  al.  1963).  Calculations  of 
A^2  are  shown  in  Figure  3  for  the  grid  described  above. 

3.0  Conclusions 

Following  Adonina  and  Shestopalov  (1963),  we  have  extended  the  set 
of  equations  for  the  reflection  and  transmission  coefficients  of  a  grid 
made  of  metal  strips  on  a  dielectric  slab.  Orthogonal  planes  of 
polarization  have  been  considered.  However,  the  equations  are 
restricted  to  the  case  where  the  plane  of  incidence  is  normal  to  the 
strips,  and  the  radiation  is  incident  on  the  side  with  the  metal  strips. 
Whereas  absorption  in  the  dielectric  substrate  has  been  treated 
rigorously  here,  absorption  by  the  metal  strips  has  only  been  estimated 
to  order  of  magnitude  following  Ulrich  et  al.  (1963).  Throughout  this 
work,  the  approximation  has  been  made  that  the  wavelength  X  is  much 
longer  than  the  grid  period  g.  The  limit  to  the  validity  of  this 
assumption  appears  to  be  ‘  >  5g  for  free  standing  strips,  but  it  is 

not  as  well-known  when  the  grid  is  on  a  dielectric  substrate.  The 
remaining  problems  will  be  addressed  in  further  work,  and  an 
experimental  evaluation  of  grid  performance  will  be  made. 
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APPENDIX  C 

Measured  Transmission  and  Reflection  of 
Grid  Polarizers  on  a  Dielectric  Film 

by 

O.A.  Simpson 

R.A.  Bohlander 

R.E.  Forsythe 

In  this  section  of  the  report  we  present  preliminary  results  of 
transmission  measurements  on  2  grid  polarizers  and  a  clear  Kapton  sub¬ 
strate  with  both  a  Fourier  Spectrometer  in  the  spectral  region  10  to 
100  cm"^  and  a  far-infrared  (FIR)  optically  pumped  laser  at  four  frequen¬ 
cies  between  8.2  and  61.3  cm~\  The  two  separate  sets  of  data  are  in 
excellent  agreement. 

The  basic  experimental  apparatus  for  the  Fourier  measurements,  a 
Grubb  Parsons  cube  interferometer  is  shown  in  Figure  C.l.  An  evacuated 
chamber  with  a  grid  polarizer  and  mirror  arrangement  (not  shown)  directed 
the  linearly  polarized  broad  band  FIR  radiation  to  normal  incidence  with 
the  second  polarizer  (or  clear  substrate).  The  grid  orientation  was  normal 
to  the  plane  of  polarization  of  the  radiation  for  maximum  signal  transmission. 
The  transmitted  radiation  then  entered  a  Golay  cell  detector.  Transmission 
spectra  were  determined  from  the  ratio  of  spectra  obtained  for  the  second 
polarizer  (or  Kapton  substrate)  and  spectra  with  no  sample.  The  measured 
values  are  shown  graphically  in  Figures  C.2-C.4. 

The  grids  tested  were  fabricated  by  Dr.  G.  Lamb  at  NASA  Goddard 

o 

Space  Flight  Center.  He  vacuum  deposited  5000  A  of  aluminum  on  Kapton  13 
Min  thick  and  then  photo-etched  the  grid  pattern  of  metal  strips.  The 
sample  of  Kapton  film  was  processed  in  the  same  way,  but  all  of  the  aluminum 
was  etched  off. 

A  block  diagram  of  the  optically  pumped  laser  is  shown  in  Figure  C.5. 

A  20  watt  CO^  laser  pumps  a  FIR  laser  of  the  metal  waveguide  type.  The 
FIR  beam  was  divided  by  the  beam  splitter  BS2.  The  reference  signal  1^, 
after  detection  and  amplification  by  a  Golay  cell  and  lock-in  amplifier, 
entered  the  denominator  channel  of  the  ratiometer  RA.  The  beam 
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then  entered  a  Fabry-Perot  interferometer  (FP)  which  was  used  both  to 
measure  the  wavelength  and  as  a  filter  to  eliminate  unwanted  wavelengths 
that  may  have  lased  simultaneously  with  the  desired  signal.  Since  the 
radiation  from  this  type  laser  is  generally  unpolarized,  the  signal  la 
was  passed  through  a  polarizer  to  produce  the  desired  plane  of  polarization. 
The  signal  was  then  transmitted  through  (or  reflected  from)  the  sample 
polarizer  and  detected  by  another  Golay  cell.  This  signal,  after  lock-in 
amplification,  became  the  numerator  signal  in  the  ratiometer.  The  ratioing 
process  was  important  in  eliminating  source  fluctuations  in  the  FIR  power 
level.  We  will  call  values  of  Ia/Ib  normalized  signals. 

The  transmission  of  a  sample  polarizer  was  determined  from  the  ratio 
of  the  normalized  signal  with  and  without  the  sample  in  the  laser  beam. 

The  reflection  coefficient  was  determined  by  the  ratio  of  the  normalized 
signal  when  the  radiation  was  reflected  from  the  polarizer  and  that  for 
reflection  from  a  mirror. 

Results  of  the  laser  measurements,  together  with  block  diagrams  of 
each  experimental  arrangement  are  shown  in  sections  I,II,III. 

CONCLUSION 

In  evaluating  the  performance  of  the  grids  we  have  used  the  theoretical 
framework  in  Appendix  B.  The  polarizers  we  have  tested  were  originally 
intended  for  millimeter  wavelengths.  However;  theory  suggested  that  they 
should  be  of  use  to  much  shorter  wavelengths  and  so  we  have  tested  them 
to  wavelengths  as  small  as  100  ^m.  It  is  apparent  from  Figures  C.3  and  C.4. 
Table  2  that  the  polarizers  perform  well  at  millimeter  wavelengths,  but 
a  close  comparison  with  theory  cannot  be  made  in  the  absence  of  good 
absorption  coefficient  and  refractive  index  data  for  the  Kapton  substrate. 

At  short  wavelengths  the  polarizer  transmission  and  reflection  are  very 
different  from  expectation.  When  the  plane  of  polarization  incident  on 
polarizer  #2  is  perpendicular  to  the  grids,  it  is  expected  that  the  trans¬ 
mission  would  be  close  to  that  for  the  Kapton  film  alone.  Instead  the 
transmission  is  about  a  factor  of  0.5  lower.  From  a  comparison  of  tables 
I  and  II  one  can  see  that  the  transmission  and  reflection  are  almost 


complementory,  so  that  the  discrepancy  does  not  arise  in  the  main  from 
absorption.  Further  experimental  and  theoretical  work  is  needed  to  un¬ 
cover  the  reason  for  the  discrepancies.  The  one  case  where  the 
discrepancies  are  not  serious  is  when  the  grid  period  is  very  small, 
as  is  the  case  for  the  data  shown  in  Figure  C.4.  The  predicted  spectrum 
for  this  case  is  very  close  to  that  of  the  substrate  and  agrees  well 
with  the  data.  The  conclusion  drawn  from  this  is  that  polarizers  for 
use  with  interferometers  should  have  very  small  grid  periods,  and  the 
period  of  10  pm  chosen  for  subsequent  fabrication  was  about  as  small 
as  could  be  fabricated  accurately  with  the  equipment  available. 

The  primary  parameters  of  the  gridded  beamsplitter  were  aperture, 

95mm;  Mylar  substrate  thickness,  6  pm;  strip  width,  approximately  4  ~m; 
and  strip  period,  10  pm.  Four  types  of  techniques  were  used  to  measure 
the  transmission  properties  of  these  grids,  and  one  to  measure  the 
reflection  properties.  Figure  C.6  shows  the  layout  of  instrumentation 
used  to  make  transmission  measurements  with  the  optically  pumped  laser 
at  frequencies  of  .525,  1.84  and  2.54  THz.  As  shown  in  Figure  C.7, 
the  transmission  of  the  polarizer  is  high  as  expected  when  the  plane  of 
the  polarization  is  perpendicular  to  the  strips.  Comparison  is  made 
with  theoretical  calculations  done  according  to  the  theory  of  Adonina 
and  Shestopalov  and  good  agreement  is  found.  On  the  basis  of  the 
work  described  above,  this  is  to  be  expected  only  when  the  wavelength 
is  much  greater  than  the  grid  spacing,  as  is  true  here.  In  this  case, 
the  polarizer's  transmission  is  mainly  limited  by  that  of  the  substrate 
as  also  shown  in  Figure  C.7.  The  transmission  of  the  opposite  polarization 
is  small  as  expected  (see  Figure  C-7),  although  not  quite  as  small  as 
predicted  by  the  theory  used.  Table  C.4  gives  values  of  the  polarization 
ratio;  i.e.,  the  ratio  of  the  signal  for  polarization  perpendicular  to 
the  strip  and  that  parallel. 
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II.  Traosaiaalon  coefficient  measurements. 


Subatrate  or 

Quart*  Fabry-Perot  p0l.  l  Pol.  2 


I.  Reflection  coefficient  measurements. 


Iris 


i  \  Pol.  2 

* 


Note:  the  translation  slide  mount  allowed  the  mirror  and  polarizer  to 

he  interchanged  ...  coincidence  of  the  two  surfaces  was  constantly 
checked  by  the  reflection  of  a  HeNe  laser.  The  second  iris  has  a 
larger  aperture  than  the  first  and  absorbing  sheets  (Eccosorb)  on 
on  both  sides  to  stop  stray  radiation. 

Polarizer  #1  has  a  strip  spacing  of  50.8  pm  and  a  strip  width  of 
12.7  urn  on  a  Kapton  substrate  12.7  uni  thick. 

Polarizer  #2  has  a  strip  spacing  of  50.8  um  and  a  strip  width  of 
19.1  um  on  a  Kapton  substrate  12.7  um  thick. 


Table  1.  Polarizer  reflectivity  data. 


X  (um) 

Line 
Pol.  1 

Orientation 

Pol.  2 

Measured 

Reflection 

570.5 

V 

V 

0.137 

570.5 

V 

H 

0.930 

185.4 

V 

V 

0.458 

185.4 

V 

H 

0.970 

1217 

V 

V 

0.045 

1217 

V 

H 

0.979 

V  -  vertical  (  X  to  page) 

H  -  horizontal  (  ll  to  page) 


FREQUENCY  (cm’1) 


Fig.  C.4  Fourier  transmission  spectra  of  a  fine  Kapton  -  aluminum  grid 
polarizer.  Polarizer  has  a  strip  spacing  of  3  yra  and  a  strip  width  of 
1.5  ym  on  a  Kapton  substrate  12.7  ym  thick. 


FREQUENCY  (cnfA) 


Fig.  C.3  Fourier  transmission  spectra  of  the  Kapton  -  aluminum  grid 
polarizer.  Polarizer  has  a  strip  spacing  of  50.8  Mm  and  a  strip 
width  of  19.1  urn  on  a  Kapton  substrate  12.7  um  thick. 


Figure  C.6  Experimental  arrangement  for  measuring  the  transmission 

of  a  mylar-aluminum  polarizer,  Pol.  2  with  an  optically  pumped  laser. 
Components  are:  B.S.,  beam  splitter;  F.P.,  Fabry-Perot  interferometer; 
Det. ,  golay  cell  detectors.  Polarizer  #1  has  a  strip  spacing  of  50.8 
microns  and  a  strip  width  of  12.7  microns  on  a  Kapton  substrate  12.7 
microns  thick. 
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Figure  C.9  Transmission  and  reflection  test  setup  using  220  GHz  Radiometer. 


A  Fourier  transform  spectrometer  in  an  arrangement  shown  in  Figure 
C.8  was  also  used  to  make  measurements  of  the  transmission  for  the  case 
of  the  polarization  perpendicular  to  the  strips  in  the  same  frequency 
range  as  Figure  C.7.  These  results  were  in  good  agreement  with  the 
laser  measurements,  though  of  somewhat  less  accuracy. 

A  radiometer  at  220  GHz  having  a  sensitivity  AT  .  of  <  0.1 °C  was 
used  to  measure  the  loss  introduced  by  the  polarizing  beamsplitter  when 
inserted  in  front  of  the  radiometer's  horn  antenna,  as  shown  in  Figure 
C.9a.  The  polarization  of  the  horn  was  perpendicular  to  the  grid 
strips.  A  loss  value  of  1  %  was  found  and  is  consistent  with  expectations 
based  on  the  laser  measurements.  Similarly,  reflection  loss  was  measured, 
for  the  opposite  relationship  between  the  polarization  and  the  strip 
orientation  as  shown  in  Figure  C-9b,  and  found  to  be  1%.  This  is 
consistent  with  expected  losses  due  to  the  transmission  of  this  polarization 
(as  may  be  estimated  by  extrapolation  from  the  laser  measurements)  and 
due  to  absorption  in  the  substrate  and  in  the  grid. 

Finally,  measurements  have  been  reported  elsewhere  of  the  polarization 
ratio  at  96  GHz.  A  value  of  630  was  found,  which  as  expected  is  a  little 
higher  than  that  found  at  220  GHz. 
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TABLE  C.4  Polarization  Ratios  of  the  Grid  Polarizer  Measured 

with  a  Laser  Source 
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Figure  C.7  Transmission  Spectrum  of  a  Mylar-Aluminum  Strip  Polarizer.  Strip 

period  10  microns,  strip  width  4  microns,  substrate  thickness  6.35  microns. 
Measurements  with  a  pumped  laser:  polarizer  strips  perpendicular  to  input 
polarization  •  ;  polarizer  strips  parallel  to  input  polarization  O  ;  Mylar 
substrate  only  A.  The  theoretical  curve  —  •  —  •  —  is  for  the  case  where 
the  polarizer  strips  are  perpendicular  to  the  input  polarization. 
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The  effects  of  the  atmosphere  have  been  paradoxically 
both  a  hindrance  and  a  boon  to  the  development  of  near- 
millimoter-wave  technology  for  manv  years,  this  spectral 
tango  was  handicapped,  not  only  bv  formidable  develop¬ 
ment  problems.  but  also  bv  concerns  about  atmospheric 
attenuation  which  is  significantly  higher  than  at  longer 
wavelengths  A  resurgence  of  interest  in  near-millimeter 
waves  <  amp  about  in  response  to  problems  at  shorter  wave¬ 
lengths  Infrared  svstPms.  while  able  to  operate  in  both  dav 
and  night,  are  frequently  unable  to  perform  through  clouds, 
tug  •  >r  -moke  Near-millimeter  waves  ire  able  >u  oenetrate 
w-se  utisr urantx  .viih  eason.mle  imionimi-i-  n  ;ang> 
and  angular  resolution  Activity  n  near-mniimetor-vvavi' 
propagation  research  has  increased  *n  recent  .ears  n  re¬ 
sponse  to  needs  tor  better  detmitions  of  the  limitations  In 
addition,  it  is  recognized  that  certain  processes  in  molecular 
phvsics  or  atmospheric  optics  can  be  studied  with  special 
Ktvantage  in  this  part  of  the  elec  tromagnetic  spectrum 

A  number  ot  <  ommencfable  'eyievvs  have  ippoared  (eg 
!l-l”i)  vvnirrt  examine  m  sortie  detail  the  nailable  literu- 
•ure  The  present  review  will  not  attempt  to  repeat  the,r 
>  r  1 1 1  c  «.i I  nitons  but  rather  provide  a  relatively  terse  account 
summarizing  what  is  known  and  what  are  the  frontier 
issues  The  review  will  cover  effects  of  importance  in  the 
frequency  range  90  to  1000  GHz  and  will  be  divided  in  two 
parts  to  be  published  separately  In  this  first  part,  the  effects 
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of  molecular  absorption  emission  and  n  trac  tion  and  -■ 
due  to  turbulence1  will  tie  covered  !  best  at.  n.  ifn.ii;v 
thought  of  as  clear-air  phenomena  hut  -is., 

important  component  of  adverse  weatne-  wtv,  • 

will  be  reviewed  in  a  later  patter 


II  MOlti  i  i  -VK  ABsc.’RBTION 

A  Thpori 

The  ma|or  molecular  absorbers  m  the  rn  ar-millimet.  i 
wavelength  region  are  HO  and  O  A  t-  v\  minor  cue 
stituents  can  be  seen  weaklv.  such  as  O  and  have  beep 
the  subiect  both  of  ground-based  and  high-  iitifud.  nv 
gations  (eg..  [8],  [9])  However  these  have  a  negligible 
impact  on  utilizations  of  the  band  tor  communications 
radar,  and  the  like  Molecular  absorption  at  tin's.-  trequen 
cies  occurs  principally  through  the-exc itation  ot  nitutmn.ii 
•runsitions  and  .*s  strength  s  isuuilv  '  ■  :  ' 

•ize  or  ’tie  molecules  >ip<  :-n  :-oo.--  ~mn 

oxygen  is  the  weaiser  ot  the  -wo  principal  it  -on  ■■ 
cause  it  is  a  homonudear  diatomic  with  no  "»  ni'on 
moment.  It  is  able  to  make  relatively  vvojk  transitions  s,nc" 
it  exists  in  a  triplet-sigma  ground  state  with  two  uncoupled 
(electron)  spins,  which  give  the  molecule1  a  magnetic  mo¬ 
ment  Absorption  due  to  bxvgen  snows  up  .  leariv  not  non1 
its  intrinsic  stremgth  Put  Tom  me1  Large  .  one  , miration  , . •  - 
-ri  (lie  atmosphere  Hie  water  moiec  uie  'tie  'tiler  •- (tv; 
has  a  strong  elector  dipole  ;nrl  m  sp.te  ot  me  gnt  itoms 
ui  the  moiec  ule  the  spacing  ot  the  .iPsorption  te'onanc  i->-. 
relatively  close  oue  to  'he  molecule  '  bent  shape1  met  the1 
consequent  asymmetry  of  its  rotational  inertia  Absorption 
resonances  of  both  oxvgen  and  water  molecules  are  shown 
in  the  calculated  spectrum  in  Fig  1  The  level  ot  absorption 
between  these  resonanc  es  has  a  .  onsideratlli1  •  ontnpution 
'roni  the  wings  ot  strong  osnnani  ns  it  v.iiec  'lignm 
trequcTic  les 

An  important  tool  in  modeling  atmospheric  aPsorptmn  is 
the  compilation  of  atmospheric  tbsorption  esonaric 
parameters  which  is  maintained  pv  the  Air  Force  Cj.'o- 
phvsics  Laboratory  ( 10],  parallel  efforts  should  also  be  recog- 
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Fig.  t.  Atmospheric  absorption  coefficient  (or  a  (lori/ontjl 
path  near  sea  level  calculated  with  (hr  model  in  jj?)  Solid 
lino  is  (ho  (oral  prodictod  for  absorption  b\  water  and 
oxvgen  molts  ules  m  tho  atmosphere  Tho  positions  ot  oxvgen 
linos  aro  indnatod  with  arrows  Tho  dashod  lino  shows  tho 
lontmuum  r  omponont  of  tho  model  laltulatod  with  (_') 
Conditions  assumed  T  =  288  K  dtv  air  prossuro  Pt  =  101  3 
kPj  ( 1013  mbar).  (tarlial  prossuro  of  water  vapor  o  =  I  '  WPa 
( I  "  mbar.  i  o  saturalod)  Attor  (3’) 


mzod  [I  1],  [12]  A  word  about  nomenclature  is  appropriate 
here  The  term  resonance,  to  denote  regions  of  peak  molor  - 
ular  absorption,  is  phvsuallv  more  satisfying,  however,  the 
convention  in  spectroscopy  is  to  call  these  peaks  "fines." 
and  this  paper  will  comply  to  avoid  contusion  with  the 
literature  The  AFGl  tables  compose  extensive  lists  ot  line 
positions  intensities,  widths,  ami  energy  levels  from  whnh 
the  absorption  loeihcient  at  a  gnen  frequency  e  may  be 
calculated  according  to 

<*(*')  =  (tLs./V.  VY,,)  O 

1  I 

where  vt  is  the  line  position.  5,(  is  the  intensity,  /(n.  e,.  yl() 
is  a  function  which  describes  the  shape  of  fhe  lines,  y,,  is 
the  linewidth.  and  q  is  the  concentration  ot  the  absorber 
As  shown,  a  simple  summation  is  made  over  all  lines  in  a 
band,  since  line-coupling  effects  mav  be  neglected  for 
atmospheric  molecules  in  the  near-millimeter  spectrum  [13] 
Allowance  is  made  for  the  temperature  dependence  ot  the 
hue  intensities  cviih  rhe  help  ot  tabulated  energy  levels.  and 
■or  -he  temperature  dependence  ot  imewidths  with  the 
help  oi  detailed  calculations  for  at  least  two  temperatures 
[14],  [IS]  Since'  the  width  ot  lines  m  the  troposphere  arises 
from  collisional  broadening,  this  parameter  depends  on  the 
barometric  pressure  and  the  partial  pressure  ot  the  absorb¬ 
ing  gas  [14],  [15]  As  a  result,  in  the  gaps  between  line's, 
absorption  bv  waler  vapor  has  twee  terms  one  that  increase's 
lieiearlv  with  water  molecule  i  onceotration  (line  broad¬ 
ening  bv  jir),  and  one  that  increases  quadratic  ally  with 
water  molee  ule  e oncentration  iselt-broadenmgi  lhe>  accu- 
rac  v  of  tabulated  data  on  line  parameters  is  reviewed  in  [  1 0] 
and  is  generally  adequate  for  modeling.  Widths  are  difficult 
to  measure,  but  by  far  the  most  uncertain  item  in  any 
model  is  the  shape  of  the  lines  a  few  linewidths  awav  from 
line  center,  particularly  when  an  entirely  theoretic al  dp- 
sc  option  is  used 

interest  in  the  physics  ot  lme>  shape  has  been  spurred  in 
the  near-millimeter  spectrum  since  much  ot  this  regions 
utility  lies  in  the  intervals  of  frequency  between  line's,  and 
there  signitie  ant  limitations  jre  imposed  by  molee  ulae  ab- 
soepbon.  even  more  so  than  in  those  regions  commonly 


used  in  the  infrared  For  many  sears  collisional  line  shape 
formulas  have  been  used  whic  h  were  derived  [  li.J  [ 1 3|  with 
a  simplification  known  as  the  impact  approximation  m 
which  it  is  assume'd  that  collisions  occur  m  mtinite-im.il 
intervals  of  time  This  c  an  give  a  tairls  ac  i  urate  pu  tun  ot 
the'  line  shape  within  intervals  of  frequenc  v  e’qual  to  a  ti  w 
times  the  linewidth  since  this  regime  corresponds  to  tune 
intervals  of  the  order  of  the  time  between  collisions  oi 
longer  (One  can  think  of  collisional  line1  broadening  i-  a 
kind  of  modulation  process)  The  details  ot  the  events 
within  collisions  begin  to  have  an  effect  on  the  line  shape 
as  the  frequency  difference  from  a  given  line  center  ap 
proaches  the  reciprocal  of  the  collision  durations  Recent 
work  has  made  progress  in  methods  to  include  the  (mite 
duration  of  collisions  in  the  theory  [19],  [20]  and  some 
general  constraints  on  the  shape'  of  lines  far  from  line 
center  have  been  clarified  For  frequenc  ie's  distant  from  a 
given  line  by  more  than  fhe  reciprocal  of  typical  collision 
durations,  some  type  of  exponential  decrease  is  required  in 
the  line  shape  formula  [21]  This  represents  lews  absorption 
tn  the  far  wings  of  lines  than  has  been  heretofore  picdu  ted 
bv  impact  approximation  formulas  and  is  understandable  as 
a  smoothing  of  fhe  modulation  process  Howeve  r  the  in 
elusion  of  a  mean  lollision  duration  is  not  the'  last  hurdle  lo 
the  development  of  an  adequate'  theory  as  mav  be  seen 
irom  an  appreciation  of  the  phvsus  ot  the  collision  proc  es- 
as  well  as  from  comparisons  of  observation  with  peedu  non 
The  problem)  faced  in  understanding  atmosphere  -pe.  t  , 
is  that  the  water  molecule  the  prime'  absorber  pi.t.-in- h- 
molecular  interae  tions  complicated  bv  relatively  -imttg  it 
tractive'  forces  which  are  to  some  degree  not  spherically 
symmetric  due  to  the  shape  of  the  molecule  Attractive 
interactions  are  believed  to  be  responsible'  tor  the  observa¬ 
tion  of  increased  absorption  between  lines  when  temper¬ 
ature  decreases  There  has  been  a  controversy  about  wheiher 
this  absorption  signals  the  formation  of  dime'rs  or  should  tic- 
understood  as  monomer  line  wings  sensitive*  to  the'  altr.u  - 
teve  part  of  the  potential  energy  ot  interaction  Fot  simpler 
molecules  (which  have  spherically  symmetric  shapes)  ab¬ 
sorption  caused  by  unbound,  metastable,  and  bound  pairs 
can  be  predicted  m  the  correct  proportions  at  a  given 
temperature  (22),  ( 23]  In  general  when  the  i.-mm  t  nu  .  - 

low  t)v  c  omparison  with  the  binding  "nergv  oi  moim  .i.ii 

jrairs  dimers  are  numerous  and  plav  an  important  . .  m 

the  absorption  ot  such  simple  molecules  Although  * n is 
condition  is  met  tor  pairs  of  water  molecules  (the  heat  < > t 
formation  is  believed  to  be  at  least  as  large  as  0  12  >A  as 
compared  with  ksT  =  0025  eV  [24]),  the  theory  ot  the 
relative  importance  of  different  kinds  ot  pair  inter, k  tions 
has  not  been  extended  as  vet  to  asvmmetrir  molee  uies 
Microwave  transitions  ot  water  dimers  have  been  ob¬ 
served  in  molecular  beams  (25 ]  and  arc-  consistent  with  the 
lowest  energy  structure  calculated  trom  molecular  orbital 
theory  (eg.,  [2b]>.  However,  molecular  beams  provide  what 
is  an  essentially  very  cold  molecular  environment,  and  an 
analysis  ot  observed  room-temperature  laboratory  spectra 
shows  no  band  structures  simply  related  to  the-  known 
low-temperature  form  of  the  dimer  |2  l)  if  dimers  in  ,  on- 
tribute  a  signific  ant  oart  of  the-  absorption  at  mom  >•  min" 
.iture  d  would  appear  that  thev  must  "vi-t  in  man',  aim 
lure's  different  from  the  lowest  energs  torn)  While  theory 
does  not  vet  predict  how  the  spectrum  ot  the  .inner  ■  vniv>  s 
with  increasing  temperature,  one  woulci  expei  t  consider¬ 
able  changes  since  current  calculations  (2b|  indicate  that 
the1  nonspherical  part  of  the  potential  energy  lunc  turn. 
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which  defines  the  particular  configuration  of  lowest  energy, 
is  relatively  weak. 

In  the  complementary  arena  of  line  broadening  theory, 
some  progress  is  being  made  toward  the  inclusion  of  more 
aspects  of  the  intermolecular  potential  energy,  but  a  com¬ 
plete  inclusion  of  what  is  known  [26]  will  be  a  formidable 
undertaking  Thus  a  satisfactory,  completely  theoretical  de¬ 
scription  of  water  vapor  absorption  is  not  presently  avail¬ 
able,  and  current  models  of  atmospheric  transmission  are 
empirical 

B  Observations 

There  ■  ;  several  fairly  recent  reviews  of  data  on  the 
transmission  of  clear  air  or  of  atmospheric  constituents 
measured  in  the  laboratory  [2],  (6],  [24],  [27],  Such  data  are 
not  easily  obtained  as  may  be  elaborated  briefly  by  the 
following  generalizations  concerning  the  various  methods 
used 

1)  Spectra  can  be  obtained  with  Fourier  spectroscopy,  but 
available  black-body,  wide-band  sources  are  weak,  and 
require  liquid-helium-cooled  detectors  and  relatively  long 
integration  times  for  satisfactory  results  Spectroscopy  per¬ 
formed  instead  with  multiple  or  tunable  narrow-band 
sources  sometimes  requires  tedious  tuning  operations,  but 
generally  offers  greater  accuracy  at  a  given  frequency. 

2)  Laboratory  studies  have  the  advantage  of  better  control 
of  the  subiect  constituents  their  pressures  and  temper¬ 
atures.  but  the  disadvantage  of  limited  path  lengths  relative 
to  the  kilometer  size  links  important  in  applications  of  the 
band  These  frontiers  have  been  pushed  back  through  ex¬ 
treme  care  in  fabrication  of  open-resonator  cells  [28]  and 
through  the  development  of  large,  untuned  resonator  cells 

129]. 

3)  In  direct  atmospheric  research  there  is  relatively  greater 
difficulty  in  specifying  the  atmospheric  conditions  and,  of 
course,  a  lack  of  control  of  these  conditions,  but  in  several 
on-going  researches,  a  greater  effort  on  characterization  is 
being  made  There  are  also  difficulties  in  obtaining  an 
absolute  scale  of  transmission  or  absorption,  since  the 
atmosphere  cannot  be  "pumped  out"  in  order  to  Obtain 
signals  from  the  same  experimental  apparatus  with  no 
atmospheric  absorption  and  thereby  to  make  allowance  tor 
instrumental  losses  and  losses  due  to  beam  divergence 
Narrow-band  measurements  have  something  of  an  ad¬ 
vantage  m  that  the  beam  divergence  in  the  far  field  can  be 
calculated  and  experiments  involving  multiple  paths  can  be 
analyzed  in  such  a  way  that  instrumental  losses  are 
eliminated  [30]  Those  who  employ  Fourier  spectroscopy, 
and  m.inv  who  use  narrow-band  sources,  must  settle  tor 
relative  transmission  measurements  from  which  the  t  hange 
ot  attenuation  tor  a  given  change  in  humidity  or  tempera¬ 
ture  can  be  obtained  When  the  database  is  large  enough, 
extrapolation  to  zero  water  vapor  is  possible  [31],  and  an 
absolute  scale  for  the  important  water  vapor  component 
can  be  established 

4)  Slant-path  radiometrv  has  also  been  used  to  measure 
emission  and  absorption,  but  the  difficulty  of  knowing  the 
atmospheric  constituent  concentrations  and  temperatures 
over  the  path  complicates  the  interpretation  of  the  results 

Despite  the  diversity  of  methods,  a  consensus  is  forming 
on  data  of  sufficient  quality  and  <  onsistencv  [2],  [6|.  [2-1]  [2’] 
that  an  adequate  basis  exists  to  construct  empirical  models 
useful  over  a  fair  range  of  clear  air  conditions  in  the  atmo¬ 
sphere.  The  spectrum  in  Fig  1  illustrates  such  a  calculation 


C  Anomalous  Attenuation  and  Recent  Work 

Several  significant  researches  have  been  reported  since 
the  reviews  mentioned  above  which  have  been  motivated 
by  another  controversy  involving  claims  for  the  existence  of 
anomalous  attenuation  ascribed  to  complexes  of  water  of 
various  sizes  [32]— [34]  The  conditions  under  which  this  is 
said  to  be  observable  are  ones  of  high  humiditv  with  clear 
air  or  fog  present,  and  ones  of  low  temperatures  in  which 
relatively  weakly  bound  complexes  might  plausiblv  be 
favored.  Consequently,  such  conditions  have  been  given 
more  attention  in  recent  work  Examples  of  anomalous 
absorption  indications  are  given  bv  a  narrow-band  labora¬ 
tory  study  in  a  special  untuned  cavity  at  frequencies  of 
115-126  and  213  GHz  [29],  [35]  and  a  Fourier  spectroscopic 
investigation  of  atmospheric  transmission  in  the  range 
150-870  GHz  [36]  Both  showed  increases  in  specific  at¬ 
tenuation  with  decreasing  temperature  (below  290  K)  which 
were  much  larger  than  expected  both  bv  comparison  with 
previously  observed  temperature  dependences  in  higher 
temperature  ranges  and  by  consideration  of  the  tempera¬ 
ture  dependences  likely  to  be  associated  with  energy  levels 
in  the  water  molecule  and  expected  intermolecular  interac¬ 
tion  energies  Moreover,  the  absolute  attenuations  derived 
from  these  experiments  were  considerable  in  engineering 
terms;  namelv,  in  the  range  5-20  dB/km  at  various  frequen 
cies  between  about  210  and  300  GHz.  values  which  it 
correct  would  have  a  noticeable  impact  on  applications  ot 
near-millimeter  waves.  A  number  of  results  nave  appeareo 
since  which  do  not  show  these  effects  ana  raise  doubts 
concerning  their  validity. 

A  careful  laboratory  study  has  recently  appeared  of  at¬ 
tenuation  at  138  GHz  [28]  An  open  resonator  was  used  m 
which  care  was  taken  to  avoid  adsorbed  lavers  of  water  on 
the  reflecting  surfaces  which  might  add  spurious  attenua¬ 
tion.  Results  were  obtained  for  temperatures  of  282  and  300 
K  and  for  various  mixtures  of  nitrogen  and  water  vapor  This 
has  been  used  to  determine  empirical  correc  tions  to  model 
spectra  (see  below)  which  have  been  used  to  make  com¬ 
parisons  with  other  data  No  support  for  low -temperature 
or  high-humiditv  anomalies  was  found  and  agreement  with 
■  fata  m  previously  mentioned  rev  iews  T 2 1  !•■>!  12  H  - 

good.  Reasonable  agreement  has  also  noon  touno  I  who 
new  results  at  1 10  GHz  [38],  [3q]  in  which  techniques  were 
used  that  were  closelv  similar  to  those  [2a1  ! ? :  1  whuh 
previously  showed  anomalous  attenuation,  although  an  al¬ 
ternate  analysis  [39]  still  finds  an  anomalous  low-tempeia- 
ture  effect  at  several  frequencies  between  3°  0  and  0 
GHz.  its  contribution  to  total  absorption  is  not  large  , 

New  Soviet  studies  have  recently  been  described  I  TO]  it'1 
of  measured  attenuation  in  the  trequencv  range  ->0-  i 
GHz  over  a  ;.5-km-long  atmospheric  path  The  zero  ot  the 
attenuation  scale  was  found  bv  extrapolation  trom  relative 
measurements  at  various  humidities,  a  method  that  was 
mentioned  earlier.  Excellent  signal-to-noise  ratios  were  ob¬ 
tained,  and  no  indication  of  anomalous  temperature  depen¬ 
dences  or  unexpected  spectral  features  were  found  m  the 
temperature  range  263  to  282  K  Previous  work  bv  ’hose 
investigators  [42]  is  consistent  with  an  absence  ot  -ui  tv 
effects,  and  recent  Soviet  laboratory  results  also  -upport 
these  conclusions  [41], 

Although  proponents  of  the  existence  of  anomalous  ab¬ 
sorption  caution  that  their  observations  may  apply  to  occa¬ 
sional  nonequilibrium  events  [32],  the  tide  appears  to  be 
running  toward  a  belief  that  such  conditions  are  either 
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quite  rare  or  nonexistent  Nevertheless,  retent  acbvitv 
spurred  by  this  controversy  has  fostered  signifit  an*  im¬ 
provements  in  millimeter-wave  propagation  models  [28], 
[37],  and  further  benefits  tan  be  anticipated 


D  Empirical  \lodols  of  Clear- Air  Transmission 


There  have  been  ma|or  efforts  to  develop  computer  mod¬ 
els  of  near-millimeter-wave  propagation  at  several  national 
laboratories,  the  U  b.  Air  Force  Ceophvsics  Laboratoiv  ]43], 
(44|.  e  g  .  the  U  S  Arms  Atmospheric  Sciences  laboratorv 
(45|,  and  the  National  Telecommunications  and  Information 
Administration  [28],  [37],  [4b]-[48]  All  have  been  extended 
to  include  effec  ts  of  adverse  weather,  which  will  be  consid¬ 
ered  in  a  later  paper 

In  the  interest  of  computing  speed,  the  general  practice 
in  these  models  is  to  compute  the  summation  in  (I)  only 
for  lines  near  the  region  of  interest  and  to  add  to  this  a  term 
that  slowlv  increases  with  frequency  This  represents  all  the 
contributions  due  to  the  wings  of  lines  at  other  frequencies 
and  due  to  other  possible  absorption  mechanisms  in  water 
vapor  which  have  been  mentioned  It  has  been  customary 
in  this  field  to  refer  tea  this  term  as  "continuum  absorption" 
although  this  is  noi  entirelv  consistent  with  previous  usage 
in  spectroscope  Its  magnitude  has  been  determined  tax 
ruling  to  observaticans  and  thus  compensates  tor  the  limita¬ 
tions  which  exist  in  current  theorv.  In  earlv  empirical  mod- 
<-lirig  efforts  [4'-a|  data  tor  limited  conditions  of  barometric 
tir.'ssuro  absolute  humiditv,  ana  temperature  were  utilized 
n  di>*  nt.  end  thus  there  was  inadequate  underpinning  to 
the  assumptions  made  concerning  the  dependence  of  the 
magnitude  of  this  term  on  those  variables.  In  particular 
inadequate  attention  was  paid  to  the  effect  of  self-broad¬ 
ening  on  the  humidity  dependence  and  to  the  steepness  of 
the  temperature  dependence  in  the  gaps  between  lines  A 
recent  tormulation  of  the  continuum  term  a,  [28],  [37] 
given  by 


=  (0  ytfePJ 


is  much  improved  in  this  regard  and  is  illustrated  in  Fig.  I 
bv  the  dashed  curve  FTere,  the  water  vapor  partial  pressure 
"  is  in  kPa  (=  '0  mbar).  the  drv  air  pressure  P  is  m  kPa. 


jM/,  mo  ■  *v* 


nnrmuum  .m 


. nt  i  s  n  iB/i.m  'n  inorrv  r  loproacn  a  .  diiiiiiu- 

im  has  been  derived  from  the  wings  or  in  empirically 
nodi  tied  line  shape  that  has  neen  successfuilv  fitted  to  data 
m  parts  of  the  millimeter-  and  near-millimeter-wave  regions 
is  well  as  to  the  infrared,  with  onlv  tour  fitting  parameters 
being  net  essarv  143], 

!  omoarisons  between  formulations  ot  c«  ,  suc  h  as  in  [.81. 
i n*  diificuif  to  made  and  assess  since  the  lines  included  m 
■  r > c ■  mi'-dv-nne  .urt  .>|  i  given  mocjel  and  the  ime  shape 
tormina  used  .ire  seldom  the  same  from  one  model  to  the 
■text.  These  choices  i.jn  have  a  dramatic  etfec.t  not  onlv  on 
the1  magnitude  or  the  complementary  continuum  term,  but 
also,  for  example,  on  its  apparent  frequency  dependence 
borne  workers  show  a  continuum  with  a  simple  frequency 
fependence  like  that  in  (2)  and  others  a  continuum  with 
'vmi  or  more  mrlec  bon  points  .'4|  [_~l  );]  The  tp- 


■pervienc  e 


trv*  *  orinnuum  is  -mii- 


,i r i *.  jttec  tec]  It.  *h**  line's  .ee-c  t*'0  lor  expiie  *1  *nc  lusion  *n  a 
modc'l  calculation  and  mav  varv  with  trequenr  v  [24]  In  the1 

>t  ontinuuio  'orrrvii.is  :  .  *  nil  to  n*mc  ’Tiber  -ri.it 

ttus  so-called  continuum  is  not  vet  a  separable  physical 


phenomenon  but  rather  an  artifice  which  is  added  m  to 
simplify  calculations  ancf  improve  agreement  witti  observa 
tions  Unnecessary  discrepant  les  with  oFiservations  (  ,m  n* 
suit  when  parts  of  a  model  arc'  obtained  from  dirt,  n  ot 
souices  and  happen  to  be  incompatible  then  i.  fl  nni.nl 
c-rable  need  for  future  standardization  ot  praetut  t ■ .  a.oi  ) 
this  source-  of  confusion  further  work  is  und.rvv.i-  u" 
comparisons  ot  the1  effectiveness  ol  c  orupc-tii'g  *  i.c r.  *  i 
modc'ls  as  rc'gards  agreement  with  otiservaiioo.  In  u.. 
course  ot  the'  present  review  it  was  non  d  th.it  tt*.  r*  .  *.* 

Sovic't  rc'sults  |40|.  [41]  mentioned  above  jitter  1mm  pi.  I  . 
tions  based  on  the  model  in  [28]  (i_|  bv  plus  ami  men.. 
dB.  km.  depending  on  the  trequenc  or  conditions  -r .  i .  •  ■: 
(sen  Note  Added  in  Proof  on  p  38) 

Besides  the  advantages  ot  computational  spend  hum  m 
striding  the  number  of  lines  that  arc'  calculated  cxplu  itlv 
improvements  have  also  been  made  m  the  speed  with 
which  regions  near  line  centers  c  an  lie  c  ale  ulated  [4  •  ]  I  In. 
has  facilitated  modeling  of  vertical  or  slant  paths  m  vvim  h 
manv  layers  are  included  in  the  c  ale  ulation  some  c  it  v\  1 1 .«  I . 
arc'  at  low  pressures  ancf  consequently  have-  v  r.  - oi r ■  ■  >w 
line's  Fig  2  is  an  c'xamplc'  ot  such  a  calculation  toe  : 
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Fig.  2.  Transmission  ot  ,i  vc'ctccat  [i.itti  tt.ci'ciit'  *i  ■ 
bianejjre!  Atmosptiere  calculated  witti  MC.I  I  > 
modc'l  1 4  < J  The'  smatf  line  .ihsncptuir  u.  .  ■  .( 
node  eahle  tietwec'ii  .*'T<  and  **c  (,Hz  I*.  ]...  *.. 

atmospheric  constituent,  .mh  ,|.  .unni-  Ml.  .  :: 
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pels  with  zenith  attenuations  v.nn  n  in 
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Emission  and  refraction  bv  the  itrio.uh. 
treated  i  .nlv  briotlv  since  orc'Oii  lion.  1  *  •  • 

..ui  tie  derived  m  urge  measure,  r.im  .  -n.'W 


losorotion  spectrum  onivioeo  *tii 


.  tiar.u  terizjtion  o!  the  mot.n'oioyv  *r  i 
is  relatc'd  to  absorption  tiv  tvin  tthoit  s  .m 
tion,  bv  the'  Kramers- Kromg  relationship  j'  I 


A.  Emission 


In  die  ni> jr-mnlime'er  uni  n  * 
■‘ierm.il  - .ii.'cgv  .mi 1 1, -. i  .  ,.n.;n!"  . 
i!  the  emitting  'utist.ii.. /•  n 
Ravlc'igh  leans  approximation  As  a 


T„,  or  antenna  tempe'atui 


LlU'MI 


P  received  from  the  atmosphere  bv  a  radiometer  at 
frequency  v  by  [51) 

p(  v)  -  LrT„(  v)v  dvAQ/c  (3) 

where  dv  is  the  receiver  bandwidth  and  4  and  (2  are  the 
effective  receiver  arpas  and  beam  solid  angles  respertneh 
If  the  atmosphere  could  be  represented  bv  one  homoge¬ 
neous  slab  of  air.  TH  would  be  given  bv 

h  -W-t)T  T{v)t  (4) 

where  t  is  the  transmission  through  the  slab.  Td  is  the  air 
temperature,  and  the  last  term  is  due  to  the  cosmic  back¬ 
ground  and  is  small 

In  (4)  it  is  assumed  that  the  atmosphere  is  in  thermody¬ 
namic  equilibrium,  and  thus  its  thermal  emissivity,  by 
Kirchhoff's  law,  is  equal  to  its  fractional  absorption  (1  -  f). 
Gradients  with  height  of  both  temperature  and  molecular 
concentration  are  modeled  by  approximating  the  integral 
relationship 

rB(»)  =  [X  TJ(z)K(i'.z)  dz  +  T,(v)t  (5) 

by  a  summation  over  thin  spherical  shells  in  the  atmosphere 
which  are  taken  to  have  uniform  conditions.  Here  K(v,z)  is 
a  weighting  function  which  determines  the  contribution  of 
the  air  temperature  at  altitude  z  to  the  sky  brightness 
temperature  and  is  given  bv 

d 

ATI  v.z)  =  —  MO. 4)  (b) 


t(O.z)  >■  10"° (7) 
if  a  is  expressed  in  dB/km. 

When  one  wishes  to  determine  the  transmission  of  the 
atmosphere  over  vertical  or  slant  paths,  one  can  in  principle 
point  a  radiometer  at  a  source  of  radiation  outside  the 
atmosphere,  such  as  the  sun.  and  correct  the  received  signal 
for  what  is  known  of  the  source  radiance  and  the  antenna 
coupling  efficiency  (both  functions  of  frequency)  In  prac¬ 
tice,  these  corrections  are  difficult  and  some  simplifications 
may  be  reaiizeb  it.  instead,  emission  trom  the  atmosphere  is 
measured  and  the  corresponding  transmission  is  interred 
[52).  If  the  distribution  of  atmospheric  temperature  with 
height  is  known,  from  radiosonde  observations  for  example, 
measured  spectra  of  Tg  can  be  inverted  through  (5)-(7)  to 
find  the  transmission  of  the  path  If  the  water  vapor  distri¬ 
bution  is  also  known  or  can  be  estimated,  a  mean  radiating 
temperature  can  be  assigned  bv  a  Curtis-Godson-tvpe  ap¬ 
proximation  [50],  [53]  Then  (4)  can  be  used  to  obtain  a 
omple  approximate  inversion  to  transmission  An  example 
of  data  inverted  bv  this  means  [52]  is  shown  in  Fig  5  and 
compared  with  a  calculation  based  on  a  recent  empirical 
model  [43]. 

This  and  other  recent  comparisons  [37],  [54],  [55]  between 
models  and  slant-path  observations  have  shown  reasonable 
agreement  The  significance  of  these  comparisons  is  tem¬ 
pered,  however,  by  the  scant  imnunt  ol  lata  available 
!  oncoming  'he  humidity  and  temperature  profiles  over  the 
paths  studied  some  workers  have  attempted  to  derive 
continuum  expressions  for  use  in  empirical  models  on  the 
oasis  ot  slant-path  emission  observations  [52]  Although 
these  do  not  greatly  differ  in  predicted  absorption  from  the 
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Fig.  1.  Slant-path  transmission  interri-d  trom  .‘mission  rr.*w- 
surements  [52]  (dots)  and  compared  with  a  model  i  jl,  ulj- 
tmn  of  transmission  [43]  Slant  seunt  =  '  1  .ntegrjmd  « 

vapor  column  densitv,  derived  trom  radiosonde  data  =  i 
mm  precipitablp  After  [43] 


models  referenced  in  the  previous  section  thev  should 
probably  be  given  less  weight  of  confidence  than  models 
based  on  laboratory  or  horizontal-path  measurements  where 
the  atmospheric  conditions  can  be  better  specified 

There  is  interest,  of  course,  in  levels  of  atmospheric 
emission  for  their  own  sake  as  thev  impart  astronomiv  at 
observations  [52],  [54],  [55],  communications  [p]  [5n]  and 
radiometric  surveillance 

B.  Refraction 

In  the  near-millimeter-wave  range  atmospheric  retractw- 
itv  can  be  described  as  a  sum  of  a  constant  term  \ 
derived  trom  a  sum  ot  contributions  bv  lines  at  ail  ireauen- 
cies.  and  a  dispersion  term  V,.  due  to  the  water  jnu 

oxygen  lines  in  the  region 

N-  N0+  N,r  (8) 

(Note:  refractivity  is  equal  to  the  refractive  index  minus  one 
and  is  usually  expressed  in  parts  per  million  of  unity,  ppm 
equivalent  to  the  term  N-units  which  is  frequently  seen  )  In 
common  with  the  radio-frequencv  and  microwave  regions 
of  the  spectrum,  the  nondispersive  term  used  in  near-milh- 
meter-wave  models  [57] 

N-  K,(PJT)-Z:'  +  K,  (  e/  T)  Z”  -  K,{  e/T)  Z„  (9) 

depends  strongly  on  the  amount  ot  water  vapor  present  n 
the  atmosphere  (expressed  here  as  the  pattiai  pressure 
This  differs  from  the  case  in  the  visible  or  near-mtrare.: 
parts  of  the  spectrum,  where  the  nondispersive  retractiv itv 
depends  almost  entirely  on  barometric  pressure  (equal  to  e 
plus  the  partial  pressure  of  dry  air  P^)  and  temperature  T 
Equation  (9)  is  not  an  exact  theoretical  torm  but.  with 
K.  -  '760  K/kPa.  K>  -  '25  K/kPa.  and  K-  -  3  '50  v  •?* 
K’  'kPa.  one  wmch  tits  [57],  [53]  observations  well  The 
factors  Z,  3nd  2,.  are  neartv  eaual  to  unitv  [33]  [59]  atirt 
correct  tor  the  nonideal  gas  relation  between  .tensitv  ino 
pressure.  It  is  interesting  that  a  small  but  statisticailv  Mgnite 
cant  discrepancy  has  been  *ound  [5']  between  the  above 
values  of  K ,  and  K j  and  those  which  best  tit  theoretical 
refractivity  calculated  from  tables  of  line  parameters  It  has 
been  suggested  [57]  that  this  indicates  either  an  error  in 
values  ot  the  strength  or  certain  kev  lines  m  the  water  vapor 
pure  rotation  bant)  or  in  -'(feu  ot  association  ::v  water 
molecules 

The  dispersive  term  is  calculated  [37],  (4bj.  [4~],  [o0],  [bij 
as  a  summation  over  lines,  similar  to  (  M.  hut  with  a  Une 
shape  derived  trom  /  Pv  application  of  the  Kramers-Kronig 
relation  A  tvpical  value  of  is  350  ppm,  and  bv  compare 
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Fig.  4.  The  dispersive  part  of  the  atmosphere  s  refractivitv 
calculated  with  the  same  model  and  conditions  as  in  Fig  1 
For  these  conditions  the  total  refractivitv  at  0  Hz  is  350  ppm 


Significant  work  has  also  been  done  in  the  microwave 
region  (eg.  [68]).  but  millimeter-wave  or  near-millimeter- 
wave  studies  have  only  recently  been  given  prominent 
attention  The  possible  propagation  effects  which  mav  he 
catalogued  are  fluctuations  m  received  intensitv  angle  ot 
arrival  propagation  delay,  frequency,  and  polarization  ot 
which  the  first  three  mav  be  expected  to  be  sigmiiianl  at 
near-millimeter  wavelengths  In  image  formation  the  phe¬ 
nomena  can  be  manifested  as  scintillation  image  dam  mg 
and  image  blooming,  of  which  the  first  two  ari  important  at 
such  wavelengths 

As  pointed  out  in  [3],  turbulence  effects  have  sometime* 
been  discounted  for  near-millimeter  waves,  since  scattering 
by  turbulent  eddies  is  a  diffractive  effect,  which  naturally  i* 
stronger  at  shorter  wavelengths  This  is  demonstrated  bv 
the  following  expression  for  the  expected  log-mtensitv  vari¬ 
ance  of  a  received  signal  [60] 


son,  the  dispersive  term  illustrated  in  Fig.  4  is  important 
onlv  in  the  vicinity  of  the  strong  lines,  such  as  those  at  557. 
752,  and  988  GHz.  Observations  of  dispersion  in  the  near- 
millimeter  range  are  thus  far  available  only  at  low  resolu¬ 
tion  from  dispersive  Fourier  transform  spectroscopy  [62], 
but  no  significant  discrepancies  with  models  have  been 
identified 

The  practical  effects  of  clear-air  refractivity  are  to  pro¬ 
duce  phenomena  which  mav  be  grouped  into  ones  result¬ 
ing  from  large-scale  structures  in  the  atmosphere,  such  as 
diurnal  and  weather  changes,  overall  temperature  and  water 
vapor  gradients  or  inversions,  tidal  effects,  and  the  like,  and 
those  involving  small-scale  structures  such  as  those  pro¬ 
duced  by  turbulence.  The  former  will  produce  varying 
propagation  delays,  beam  bending,  beam  ducting,  multi- 
path  interference,  etc.  Except  in  the  immediate  vicinity  of 
Strong  water  vapor  lines,  as  mentioned  above,  these  effects 
should  be  closely  similar  to  those  found  in  the  microwave 
region,  since  the  refractivities  in  the  two  regions  are  very 
nearly  equal  and  the  sizes  of  the  atmospheric  structures  are 
large  compared  with  the  wavelength  in  both  cases  One 
distinction  may,  however,  be  found  for  multipath-type  ef¬ 
fects  near  the  ground  where  reflected  signals  off  the  ground 
may  contribute  to  a  propagated  wave.  In  this  case,  tne 
ground  effects  will  differ  at  near-millimeter-vvave  frequen¬ 
cies  from  those  at  lower  frequencies  due  to  significant! v 
different  scattering  properties.  Little  experimental  attention 
has  thus  far  been  given  to  these  large-scale  phenomena  in 
near-millimeter-wave  propagation.  Studies  are  being  made 
of  the  effects  of  small-scale  phenomena,  however,  as  de¬ 
scribed  m  the  following  sections 

IV  TuRBUlENCt 

When  the  wind  carries  inhomogeneities  in  the  atmo¬ 
sphere  through  a  propagation  path,  fluctuations  in  several 
properties  of  the  propagation  result.  In  clear  air,  turbulence 
is  the  cause  of  mhomogeneities  in  both  the  refractive  index 
and  the  absorption  coefficient  of  the  air.  The  refractive 
variations  explain  the  familiar  visible  phenomena  ot  star 
twinkling  and  the  shimmering  of  distant  images  on  a  hori¬ 
zon  on  a  hot  day  Vtanv  aspects  of  the  effects  at  visible  and 
near-infrared  wavelengths  have  been  explored,  after  semi¬ 
nal  theoretical  work  bv  Rvtov  [63],  Chernov  [64]  and 
Tatarski  [65],  and  detailed  reviews  are  available  [5],  [66],  [67] 
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in  which  k  is  the  wavenumber  2 ir/\.  L  is  the  propagation 
path  length,  L0  is  the  size  of  the  largest  turbulence  eddies, 
and  C„  is  the  refractive  index  structure  parameter  a  mea¬ 
sure  of  the  spatial  variance  of  refractive  index  The  reason 
that  turbulence  effects  cannot  be  neglected  in  the  near-mil- 
1 1 meter- wave  picture  is  the  significant  contribution  bv  water 
vapor  to  the  refractive  index  that  distinguishes  this  regime 
from  shorter  wavelength  regions,  as  mentioned  >'jriier  and 
the  tact  that  water  vapor  mhomogeneities  n  the  atmo¬ 
sphere  can  be  large  Fluctuating  signal  strength  is  not  the 
only  propagation  issue  of  interest  Since  at  these  wave¬ 
lengths,  high-precision  tracking  and  adequate  image  forma¬ 
tion  systems  may  be  constructed  with  antennas  of  mod 
erate  size,  it  is  particularly  of  interest  to  know  whether 
fluctuations  in  angle  of  arrival  may  show  up  as  significant 
[69], 

The  near-millimeter-wavelength  regime  is  of  consider¬ 
able  interest  for  a  number  of  other  reasons  Some  con¬ 
troversy  has  arisen  over  the  interpretation  of  observations 
near  the  center  of  water  vapor  lines  of  moderate  strength 
where  the  refractive  index  must  be  treated  as  c  omolex  i* 
will  be  discussed  further  m  later  sections  Mien  'renuenc  e-* 
are  of  practical  interest  as  possible  choices  tor  .  ontroiled- 
range  communication  links.  Another  distinction  trom  shorter 
wavelengths  is  that  near-millimeter  wavelengths  are  no 
longer  much  smaller  than  turbulent  eddies,  but  rather  are 
comparable  with  the  smallest  ones  This  so-called  inner 
scale  of  turbulence  r',,  is  of  the  order  of  a  few  millimeters 
On  the  other  hand,  it  is  possible  to  construct  antennas  or 
antenna  arravs  which  are  as  large  as  the  largest  eddies,  or 
the  outer  scale  ot  turbulence  L  ,,,  which  range*  trom  a  meter 
to  a  few  meters  m  size  tor  propagation  near  the  grouno 
and  possible  to  have  propagation  ranges  tor  which  the 
Fresnel  zone  of  size  yX7  is  also  comparable  with  Lt,  This 
has  given  investigators  an  opportunity  for  interesting  tests 
and  extensions  of  theoretical  models 


4  Theory 

If  one  may  assume  that  turbulent  eddies  are  distributed 
in  the  neighborhood  of  a  propagation  path  in  a  homoge¬ 
neous  and  isotropic  fashion,  the  Kolmogorov  model  [6ts|  ot 
the  eddy  spectrum  says  that  the  refractive  index  spatial 
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variance  is  proportional  to  the  two-thirds  power  of  the 
separation  of  two  measurement  points  A r  when  A r  is 
between  fn  and  Ln,  i.e. , 


The  brackets  denote  ensemble  averages  which  can  usually 
be  replaced  by  time  averages  (according  to  an  ergodic 
hypothesis)  The  constant  of  proportionality  C'„,  which  also 
appeared  in  (10),  is  the  key  refractive  index  structure  param¬ 
eter  relating  the  turbulence  structure  to  the  resulting  opti¬ 
cal  effects.  In  general,  the  index  of  refraction  must  be 
considered  to  be  complex,  and  strictly,  terms  should  be 
added  to  (10)  for  the  contributions  of  the  variance  of  the 
real  and  imaginary  parts  and  for  the  covariance  of  these 
components.  The  terms  involving  the  imaginary  part  are 
due  to  spatial  variations  in  absorption,  and  it  has  been 
estimated  that,  except  in  special  circumstances,  the  contri¬ 
bution  of  absorption  to  fluctuations  will  be  difficult  to 
observe  in  the  near-millimeter-wavelength  range  [60],  [61], 
The  discussion  will  return  later  to  this  topic,  but  for  the 
time  being  can  be  simplified  by  ignoring  the  terms  involv¬ 
ing  the  imaginary  part  of  the  refractive  index 

Since  both  air  density  and  water  vapor  concentration 
affect  the  refractive  index,  the  structure  parameter  C;,  can 
be  expressed  in  terms  of  analogous  structure  parameters  for 
temperature,  humidity,  and  barometric  pressure  Pressure 
variations  can  generally  be  neglected  [60],  and  C„  is  given 
bv 


Ci  -  At- 


-r  2A,Aq 


<r>(Q ) 


where  the  A  coefficients  have  been  calculated  from  the 
refractive  index  spectrum  [60],  [61]  and  Q  denotes  water 
vapor  concentration  Frequently,  the  earths  surface 
evaporates  or  transpires  a  considerable  amount  of  water 
vapor  In  these  instances,  the  humidity  fluctuations  in  the 
boundary  layer  are  considerably  larger  in  relative  terms  than 
the  temperature  fluctuations  [60],  [70],  say 

a 

(Q)  (TV 

mil  humiditv  variance  dominates  the  retractive  index  struc  ¬ 
ture  This  can  result  in  values  of  the  near-miliimeter  wave 
which  are  larger  bv  an  order  of  magnitude  or  more  than 
ones  in  the  infrared,  which  are  mainly  sensitive  to  tempera¬ 
ture  structure 

That  this  is  not  always  the  case  mav  be  seen  from  the 
following  considerations  The  last  term  in  (12)  is  propor¬ 
tional  to  the  covariance  of  temperature  and  humiditv  C ,r, 
=  AFAQ)/Ar-' '  whic  h  can  generally  he  approximated  tor 
the  boundary  layer  bv  [60] 


(T)(Q) 


Cq  £v 

(QY  ( r>’ 


The  positive  sign  is  usually  found  during  the  day  and  the 
negative  sign  during  the  night,  due  to  a  change  in  the 
direction  of  the  temperature  gradient  near  the  ground 
Since  4 ,  is  negative  tie.,  retractive  index  decreases  with 
increasing  temperature),  the  last  term  in  (12)  acts  to  reduce 
the  near-millimeter-wave  refractive  index  structure  Under 
some  drv  conditions,  it  can  substantially  cancel  the  contri¬ 
bution  of  the  other  two  terms  [71]  and  result  in  a  situation 
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for  which  the  value  of  C2n  is  smaller  in  the  near-millimeter 
region  than  in  the  near  infrared 

A  formulation  was  given  earlier  in  (10)  for  the  expected 
level  of  intensity  fluctuations  that  applies  [60]  when  I)  the 
wavefronts  are  spherical,  2)  the  transmitting  and  receiving 
antennas  are  small  relative  to  a  Fresnel  zone  AT.  and  in 
turn.  3)  the  Fresnel  zone  is  small  compared  with  the  outer 
scale  Ln  The  corresponding  expression  for  angle  of  arrival 
variance  is  [64] 

•i-0.54tC£p-,/J  (14) 

where  p  may  be  interpreted  as  the  spacing  of  small  anten¬ 
nas  (or  the  width  of  a  large  antenna)  used  in  measuring  the 
angle.  Since  with  two  or  more  antennas  in  an  arrav,  one 
actually  measures  the  phase  difference  variance,  known  as 
the  phase  structure  function  D*,  an  alternate  description  is 

D„,  -  (A<f>2)  -  1.09k2iC;p5/3  (15) 

where,  in  both  (14)  and  (15),  the  dependence  on  p  holds 
when  AT  <  p  <  L„  [72].  At  larger  antenna  spacings.  one 
can  sense  the  outer  scale  by  a  drop-off  in  the  slope  of  D0 
with  p  A  transverse  coherence  length  p^  is  sometimes  used 
which  is  the  value  of  p  for  which  D0  =  2  rad  Since  in  the 
near-millimeter-wave  range  antennas  are  usuallv  smaller 
than  or  comparable  with  p^,,  the  main  source  of  imagp 
degradation  expected  from  turbulence  is  from  image  danc¬ 
ing  rather  than  image  blooming. 

In  early  derivations  of  expressions  such  as  <  1 0).  < 1  M,  and 
(15).  it  was  assumed  that  scattering  bv  eddips  was  weak 
(i.e.,  that  scattering  angles  were  small)  and  to  remtorce  this 
assumption,  it  was  additionally  assumed  that  «  f*.  or 
more  simply,  that  A  •«  f0.  Whereas  these  inequalities  are 
readily  satisfied  at  visual  wavelengths,  they  are  not  so  in  the 
near-millimeter  range,  where  A  -  (0.  It  is  reasonable  to  ask 
whether  small  eddies  are  able  to  sustain  significant  reso¬ 
nant  scattering  in  this  regime,  similar  to  that  familiar  in  Mie 
scattering  from  rain  drops.  However,  recent  work  [73]-[~5] 
has  shown  that  the  formulas  derived  under  the  assumption 
of  weak  scattering  are  still  a  good  approximation  when 
A  >  provided  that  A  <  p*,  or  L„. 

An  important  dimension  underlying  the  phenomena  dis¬ 
cussed  thus  tar  is  the  spatial  spectrum  ol  turbulent  -’ddies 
and  the  resulting  time  spectra  ot  modulations  in  propa¬ 
gation.  Between  the  inner  and  outer  scales  ot  turbulence, 
the  Kolmogorov  model  gives  [67] 

<M«)-  0.033 (16) 

for  the  spectrum  of  refractive  index  as  a  function  of  spatial 
wavenumber  k.  For  smaller  scales  (larger  wavenumbers  than 
2‘rr/^n).  refractive  index  variations  are  small  and  a  Gaussian 
tail  can  be  added  to  (16)  [65],  [67];  tor  larger  scales  (smallei 
wavenumbers  than  2ir/Lrl),  turbulent  phenomena  merge 
with  the  large-scale  variations  in  the  atmosphere  that  were 
enumerated  in  Section  lll-B  and  are  ill-defined  in  torm  In 
general,  the  variations  in  refractive  index  continue  to  in¬ 
crease  in  size  as  the  wavenumber  decreases  below  k  — 
2ir/Ln  but  with  some  decrease  in  slope  until  some  consid¬ 
erably  smaller  wavenumber  is  reached. 

When  a  transverse  wind  carries  the  turbulence  strut  turn 
through  a  propagation  path,  the  spatial  spectrum  helps  to 
determine  the  resulting  modulation  spectrum  for  propagat¬ 
ing  electromagnetic  waves.  An  approximation  that  is  fre¬ 
quently  invoked  is  that  turbulent  eddies  are  "frozen  in". 


Cwi 
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i.e.,  last  for  times  which  are  long  by  comparison  with  the 
time  taken  to  cross  the  path.  Not  all  eddies  are  equally 
effective  in  creating  modulation,  and  so  fluctuation  spectra 
measured  in  a  propagation  experiment  do  not  closely  paral¬ 
lel  (16)  For  instance,  eddies  approximately  equal  in  size  to 
the  Fresnel  zone  are  particularly  effective  in  producing 
intensity  fluctuations,  and  the  resulting  spectrum  of  such 
fluctuations  does  not  fall  off  rapidly  with  increasing 
frequency  (  until  f  >  f0  —  v(2 ir\L)~ l/J  [76],  as  illustrated  in 
Fig  5(a),  where  v  is  the  cross  wind  speed.  For  fluctuations 
in  the  phase  difference  between  two  receivers  separated  by 
p,  eddies  of  a  scale  similar  to  p  are  particularly  effective, 
and  the  fluctuation  spectrum  has  the  form  illustrated  in  Fig 
5(b).  In  contrast,  the  total  phase  in  a  link  fluctuates  slowly 
with  a  spectrum  dominated  by  the  largest  scales  in  turbu¬ 
lence  and  by  the  other  large-scale  phenomena  which  have 
been  mentioned.  For  this  reason,  it  is  not  meaningful  to 
give  an  expression  analogous  to  (15)  for  the  overall  phase 
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Fig.  5.  Theoretical  form  ot  fluctuation  soectra  lor  thp  urnii- 
igation  or  -pherirai  electromagnetic  waves  [ 7e>|  iai  l  <>g- 
■upfisitv  spp<  tral  density  normalizpd  hv  multiplying  1 IV  inp 
tat  tor  r0/orn,  Frequency  scale  normalized  bv  dividing  the 
fluctuation  frpquencv  T  bv  f0  —  v(2jtA L)~'/!  where  ;  = 
cross  wind  speed  The  dashed  line  shows  expected  low- 
frequency  behavior  when  absorption  adds  significant 
fluctuations  The  exact  frequency  where  it  |Oins  the  solid  Imp 
will  depend  on  relative  important  e  of  absorption  and  retrac  - 
turn  [41  I  (bl  Phase  difference  spec  tral  density  normalized  bv 
leveling  bv  the  phase  structure  'line  tion  D,  Frequent  \  s(  ii>- 
normalized  bv  dividing  the  fluctuation  frequency  r  bv  u 
where  p  is  thp  spacing  ot  the  points  with  different  phase  (c  i 
Total  phase  fluctuation  spec  tral  density  normalized  bv  multi¬ 
plying  bv  'he  factor  If  n;„,  The  'luctuation  frequency  scale 
is  normalized  as  in  (a) 


variance  [67],  In  summary,  for  Fresnel  zones  which  are 
typically  one  to  a  few  meters  in  size  and  for  receiver 
antennas  or  arrays  of  the  order  of  one  to  a  few  meters  m 
size,  the  important  fluctuation  frequencies  are  expected  to 
be  in  the  low  part  of  the  audio  spectrum  m  clear  air 
turbulence 

B  Observations 

The  preponderance  of  experimental  efforts  have  thus  tar 
been  directed  at  intensity  fluctuation  effects  as  these  re¬ 
quire  less  apparatus  The  frequency  dependence  [7C]  [  ] 

and  path  length  dependence  [70]  in  (10)  have  bpen  verified 
by  simultaneous  observations  at  multiple  frequencies  m 
ranges  covering  36  to  230  GHz  and  over  paths  of  different 
lengths.  Expected  changes  when  the  outer  scaled,  happens 
to  become  smaller  than  a  Fresnel  zone  have  also  been 
observed  [77],  less  attention  heretofore  has  been  given  to 
verification  of  the  refractive  index  factor  m  (10)  due  to  the 
requirement  for  a  considerable  investment  in  micro- 
meteorological  instrumentation  to  enable  predictions  ot  G 
to  be  made.  This  need  was  perceived  in  recent  years  [t  w] 
and  several  significant  efforts  are  now  underway  A  com¬ 
parison  of  observed  microwave  values  of  C:,  with  ones 
calculated  from  measured  values  of  C-Q.  C).  and  G-,  has 
shown  reasonable  agreement  [71]  Observations  have  re¬ 
cently  been  obtained  for  near-millimeter-wave  interi-m. 
fluctuations  (in  the  range  116  to  ’~3  GHz)  as  .yell  is  •,.< 
detailed  micrometeorological  parameters  ['"]  "-.].  mo  ■ 

direct  test  ot  (10)  will  thus  be  possiole 

Observations  of  the  spectrum  of  intensity  scintillation 
have  also  netted  interesting  results  The  general  form  ot  Fig 
5(a)  has  been  confirmed  [79],  including  the  expected  talloff 
with  fluctuation  frequency  to  the  -8/3  power  at  high 
frequencies  [76],  Studies  in  the  neighborhood  ]80],  [81]  ot 
the  oxygen  absorption  band  near  55  GHz,  somewhat  below 
the  near-millimeter  range,  have  demonstrated  that  mho- 
mogeneities  in  absorption  in  the  atmosphere  produce  an 
expected  increase  [81]-[83]  in  fluctuations  at  low  fluctua¬ 
tion  frequencies,  above  those  due  only  to  the  real  part  ot 
the  refractive  index,  as  sketched  bv  the  dashed  line  m  Fig 
3(a).  One  would  expect  -imilar  otter  ts  to  tie  ■•nse'v-i  :  -  * 
significant  water  vapor  absorption  lines  ,n  the  neir  m-imi.' 
ter-wave  region. 

Observations  have  been  reported,  however  ot  --tiects 
which  have  the  opposite  sense,  namely,  striking  decreases 
(bv  a  factor  of  five)  in  the  relative  intensity  fluctuation 
levels  in  the  center  of  a  moderately  strong  water  vapor  line 
near  325  GHz  [841,  [85],  as  comoarea  with  nearbv  regions 
between  lines.  Attempts  at  theoretical  treatments  i:'l  j-'t-i 
have  taken  account  ot  the  :aet  that  propagation  paths  navi- 
to  be  much  shorter  m  sue  n  instances  due  to  >n<»  'ugn  t.-v.- 
ot  absorption  (30  dB  or  more)  and  thus  retractive  y<  attenng 
can  become  relatively  more  important  than  iittractiv 
scattering  Although  it  has  been  shown  that  in  particular 
circumstances,  a  decrease  in  fluctuations  can  occur,  present 
theories  predict  much  smaller  decreases  than  those  ob¬ 
served  [85|,  [86|  A  -mail  ieoarture  ‘vis  also  been  T-onm-d 
;87j  Tom  the  usually  ••spotted  ig-uorrn.il  ii-iniiiitn-n  d 
intensitv  tor  a  region  ot  nigh  lOsornnon  n  the  milkmen-' 
wave  spectrum  There  is  an  opportunity  here  tor  turthi  r 
work  to  consolidate  the  on  turn  ot  phenomena  'h  i;  n  . .  >tv 
pany  high  absorption  which  are  cat  more  than  academic 
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interest  due  to  potential  applications  in  the  communica¬ 
tions  field,  for  example  In  further  investigations,  it  would 
be  advisable  to  eliminate  some  of  the  existing  uncertamtv 
through  observations  supported  by  detailed  micrometeorol- 
ogy  of  the  temperature  and  humiditv  structure 

Relativelv  little  attention  has  been  given  until  recentK  to 
the  important  topic  of  angle  of  arrival  or  phase  difference 
fluctuations.  Apart  from  measurements  conducted  at  micro¬ 
wave  frequencies  (68),  there  were  some  observations  near 
150  GHz  based  on  the  use  of  a  receiver  which  switches 
between  two  apertures,  which  showed  angle  of  arrival 
fluctuations  with  a  standard  deviation  of  about  '00  [trad 
[88]  The  present  authors,  with  collaborators  from  the  Na¬ 
tional  Oceanic  and  Atmospheric  Administration  and  the 
U.  S.  Army  Atmospheric  Sciences  Laboratory,  have  been 
engaged  in  a  new  study  in  the  frequency  range  116  to  173 
GHz  that  covers  both  intensity  and  angle  of  arrival  effects 
[78]  It  is  being  conducted  over  extremely  flat  and  uniform 
ground  cover  to  eliminate  as  many  complications  as  possi¬ 
ble  to  the  turbulent  structure  Simultaneous  measurements 
of  the  micrometeorological  parameters  that  characterize 
humiditv  and  temperature  structure  are  being  made  along 
with  observations  of  standard  meteorological  parameters 
Near-  and  mid-infrared  scintillations  are  also  recorded. 
Measurements  have  been  conducted  in  a  wide  range  of 
weather  conditions  throughout  one  year,  but  are  at  the 
erne  of  writing  in  oniv  a  verv  preliminary  state  of  reduction, 
foe  ingle  of  arrival  measurements  are  derived  *rom  a  10- 
•m-wiop  ov  "  I-m-men  interferometric  receiver  arrav,  and  an 
example  of  resuits  typical  of  summer  daytime  conditions 
dominated  by  humidity  fluctuations  is  given  in  Fig  6  as  a 
function  of  horizontal  receiver  spacing.  The  Kolmogorov 
5/3  power  law  (15)  is  shown  for  comparison,  and  one  can 
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see  clearly  the  expected  departure  at  larger  spacmgs  due  to 
the  effects  of  the  finite  size  of  the  outer  scale  of  turbulence 

[72J 

In  the  longer  term,  such  studies  will  be  useful  m  af¬ 
fording  the  opportunity  for  detailed  comparisons  with  the- 
orv  but  for  the  present  thev  have  also  provided  useful 
orientations  concerning  the  qualitative  importance  of  vari¬ 
ous  types  of  weather  in  producing  degradation  in  propa¬ 
gation  due  to  turbulence  (or  turbidity  trom  suspended  or 
falling  hvdrometeorites)  Summertime  conditions  when 
transpiration  or  evaporation  of  humiditv  are  strong  have 
provided  the  iargest  fluctuations  yet  seen  Over  a  path  ot 
about  1  4  km,  intensity  fluctuations  up  to  40  percent  peak- 
to-peak  and  angle  of  arrival  fluctuations  of  a  few  hundred 
microradians  have  been  observed.  At  the  other  end  of  the 
spectrum  have  been  fog  conditions  in  stationary  air.  for 
which  fluctuation  levels  have  been  about  thirtv  times  lower 
[78]  Generally,  precipitation  events  have  provided  inter¬ 
mediate  levels  of  fluctuations  This  may  at  first  seem  surpris¬ 
ing  since  during  a  rain  storm,  for  example,  there  is  plentv  of 
available  humidity  and.  in  addition  fluctuations  from 
scattering  by  rain  Offsetting  those  influences  which  in¬ 
crease  the  effect  of  turbulence  is  a  tvpical  decrease  due  to 
cloud  cover  and  the  removal  of  strong  heating  at  the 
ground.  Even  so,  of  the  factors  which  increase  the  fluctua¬ 
tions  during  a  storm,  preliminarv  indications  appear  to  rank 
the  humiditv  effects  as  being  the  most  important  ("S|  pretv. 
much  in  line  with  expectation  I3°l.  This  stresses  the  .mpor 
fance  of  the  studv  of  turbulence  not  oniv  n  ,  <ear  iir  :<u’ 
aiso  as  it  is  superimposed- on  adverse  weather 

The  discussion  of  turbulence  effects  in  this  review  has 
been  confined  essentially  to  horizontal  paths  fairlv  near  the 
ground  where  the  atmospheric  conditions  can  be  reason¬ 
ably  unifo.  m  and  readily  characterized  bv  in  situ  measure¬ 
ments.  Little  attention  has  been  given  to  turbulence  effects 
on  near-millimeter-wave  propagation  over  slant  paths  One 
might  argue  that  slant-path  effects  will  be  smaller  m  most 
regards  due  to  shorter  interactions  with  the  lasers  nearest 
the  ground  where  the  refractive  index  structure  constant  is 
expected  to  be  highest.  The  foregoing  discussion  has  also 
been  restricted  to  spherical  wavefront1-  as  •hese  ire  ofc  e 
•encountered  in  practice.  How.'vor.  nuclei-  ••  :;ur* 
and  beam-wave  rases  also  are  available  te.g.  -7]  ; 

V.  Concluding  Remarks 

An  understanding  of  the  effects  of  the  clear  itmospnev- 
on  near-millimeter-wave  propagation  is  fundamental  \ 
.vider  understanding  ot  effects  n  all  woes  )t  wi-jv  •  - 

tew  ot  the  salient  needs  'or  'urttier  .vors  .-.nun  ...--e 
discussed  m  the  preceding  sections  vui  be  outlined 

'■)  Present  models  ot  nonresonant  absorption  :  us  1 
vapor,  the  so-called  continuum  between  the  fines  ire  em¬ 
pirical.  A  breakthrough  to  an  a  priori  description  ot  the  role 
of  the  molecular  attractive  forces  of  wafer  in  collision 
broadening  of  absorption  lines  and  in  dimer  formation 
■would  be  of  maior  consequence  it  would  orovide  i  be>te- 
mderstanding  it  hvdrogen  bonding  ’>'r  i  tec  mo„><j.e  > 
veil  is  better  predictions  >r  itmospnerir  ji>s.>ro;i.»n  . 
unction  of  temperature 

2)  Even  within  empirical  aoo'oaches  to  predictions  ot 
itmosphenr  apsorotipn  them  -  -ocm 
Much  needless  confusion  arises  bv  a  lack  ot  -.tandardizat'on 
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over  the  definition  of  what  nearby  lines  to  include  in 
explicit  line-by-line  calculations.  There  is  little  on  which  to 
base  a  choice,  but  the  procedure  used  by  the  investigators 
at  the  U  S.  Air  Force  Geophysic  Laboratory  [43]  has  the 
advantage  that  it  is  useable  over  a  wide  spectral  range  An 
effort  to  bring  about  standardization  and  a  consolidation  of 
present  observations  and  predictions  into  a  format  where 
more  thorough  comparisons  can  be  drawn  would  be 
welcomed 

3)  The  discrepancy  between  theory  and  the  observed 
temperature  dependence  of  the  radio-frequency  refractivity 
is  tantalizing  [57]  The  phenomenology  of  large-scale  refrac¬ 
tive  effects  at  near-millimeter  wavelengths,  such  as  beam 
bending  and  multipath  propagation,  have  so  far  been  given 
little  attention 

4)  Greater  prominence  is  now  being  given  to  studies  of 
turbulence  effects.  Several  interesting  confirmations  of  the¬ 
ory  have  been  completed  and  many  more  comparisons  are 
in  the  offing.  One  puzzle  which  has  appeared  is  the  report 
of  a  strong  diminution  of  intensity  fluctuations  near  the 
center  of  a  line  at  325  GHz  (e  g  ,  [85]) 

Important  effects  in  adverse  weather  will  be  covered  in  a 
future  publication 

Note  Added  in  Proof 

The  author  recently  received  a  preprint  of  an  article  bv 
H  |  liebe  m  which  revisions  of  the  model  in  [28],  [37]  are 
prevented  The  improvements  are  based  on  revised  line 
data  and  a  more  refined  analvsis  of  data  obtained  in  his 
laboratory.  Thev  appear  to  bring  the  model  into  substan¬ 
tially  better  agreement  with  the  Soviet  data  [40],  [41] 
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Abstract 

The  need  for  various  co~oonents  in  millineter- 
submi  11  imeter  instrumentation  nas  led  us  to  an  in¬ 
vestigation  of  several  quasi -ontical  technioues  in 
this  wavelength  reaion.  Results  of  studies  of 
laminated  TPX  lens  material.  Eccosorb,  an  effi¬ 
cient  nuarter  wave  olate  at  225  GHz,  and  metal 
grid  oolarizing  beam  solitters  are  oresented. 

Quarter  Wave  Plate  fpr  a  225  GHz  Source 

The  efficiency  of  a  sanohire  Quarter  wave  olate 
for  a  225  GHz  source  deoenas  primarily  on  tne 
solution  of  two  general  problems: 

(1)  a  need  for  anti  -  re' 1  ec  t  ion  coatina;  the 
average  refractive  incex  of  saocnire  is 
3.24  [1]  ana  -or  the  case  o'  a  tmck 
parallel  ulane  crystal  tns  value  -orre- 
: oonos  to  a  men  sm.gie  :ur-ace  •'e'lec- 
tivitv  of  30  . 

(21  The  birefriqent  cotical  ornnertv  An  - 

n  -  n  of  the  crystal  rust  ^e  determined 
a?  225°GHz  to  establish  the  nlate  thick¬ 
ness  . 

Initial  far-infrared  laser  -easurements  at 
\=121?  uni  ( 244  GHz)  were  performed  on  a  plane 
oarallel  olate  of  sanohire  crvstal  cut  with  its 
ontic  ans  oarallel  to  tne  Plate  faces.  The 
crystal  thickness,  d  =  3.00  "h  was  such  that 
the  plate  was  close  to  a  fifth  order  ouarter  wave 
plate  at  225  GHz,  based  on  values  of  n  and  n0 
extraoolated  from  nreviouslv  determined  values  for 
frenuencies  greater  than  500  GHz  [1-3].  *he 
laser  easurements  ’nm eaten  a  significant  re¬ 
duction  n  reflection  tosses  obtained  with  *iims 
of  ''vlar  (TV,  Duoont).  n  5  1.7. 

The  difficulty  in  constructing  suen  coatings 
is  in  eliminating  the  formation  of  air  layers  or 
oockets  between  the  crvstal  substrate  and  film. 

The  air  layers  were  satisfactori  1v  removed  bv 
first  cleanino  the  surfaces  carefully  and  then 
oressinn  the  ’-'vlar  *b  "he  crvstal.  Doth  films 
were  held  in  nlace  with  neoorene  gaskets  a.nd 
damns  near  'he  outer  edge  or  the  olate.  'he 
snacino  oetween  the  Mvlxr  and  sanonire  was  men 
evacuated  ind  sealed.  With  this  method  a  maxi¬ 
mum  back  reflectivity  of  about  1.3.  was  achieved 
at  225  GHz. 

To  determine  the  value  of  in,  measurements 
were  made  of  the  effect  of  the  sanoire  plate  on 
plane  polarized  radiation,  'ne  laser  was  used 
and  its  freruenrv,  244  ghz  vas  sufficiently  close 
to  325  GHz  that  the  chanoe  in  An  between  the  two 
frequencies  was  negligible.  Significant  deviation 
front  quarter  wave  plate  conditions  were  expected 
at  the  Ijser  frequency,  'he  elliniticitv  of  the 
polarized  radiation  resuitinn  mom  two  nasses 
through  the  sanohire  plate  was  determined  and 


compared  with  theoretically  calculated  cnanf-es  in 
polarization  for  various  An.  'ne  value  \n  = 

0.345  best  fit  the  measurements  ana  is  in  n.ite 
gooa  agreement  with  extrapolations  'mm  nrs  -t  s 
results  [1-3]. 

3ased  on  this  value  of  In.  a  iawmr  r  a  v 
wave  olate  was  tonstructed  with  a  nlate  t.  i.-  ass 
of  4.33  mn.  Hgrlectivi  tv  and  ellintit’tv  '  easve- 
ments  similar  to  those  performed  with  me  -  IP 
laser  were  -ace  witn  an  Extended  Interaction 
Oscillator  at  tne  mopuenev  225  GHz.  It  ..as  fo-rc 
that  39  of  tne  •‘aaiation  after  two  transits 
through  tne  quarter  wave  nlate  had  a  -p i a za  t • 
orthogonal  to  tne  innut  ool ari za t i on . 

.ens  “ater'ai  ~or  a  215  :idip-e tzr 

in  me  course  c*  cesinninn  ano  • aprlcatinn  a 
norn-lens  antenna  'Or  a  215  In z  radiomef  r ,  it  was 
determined  tnat  aosorntion  losses  tw  the  lens 
material  Fexolite  were  substantially  higher  at  315 
GHz  tnan  at  lower  frenuencies  [4.5].  :exolite,  a 
cross  linked  nolvstvrene,  has  been  a  popular 
choice  cor  lenses  used  around  90  GHz  [6]  ••owever. 
comparison  transmi ss ion/ i nserti on  loss  .  easurements 
with  a  far-infrared  ooticallv  numoed  laser  nave 
shown  tnat  TPX  or  nol v(4-me thvl nentene- 1  gives 
much  less  absorotion  throughout  the  far-mmared 
region.  Results  of  these  laser  measurements  are 
given  in  Table  ! . 

TPX  is  clear,  slightly  amber  in  color,  and  nas 
■i  refractive  irdex  at  “illi-eter  wave  lengths  .nicn 
is  :!use  to  tnat  -o  "ne  visible  -.•••ectrvn  ■"7,2] 

A!  thougn  TPX  is  rot  as  rigid  a  “atenai  as  Fexo¬ 
lite.  'here  are  no  serious  nroolems  in  -acnimr.a. 
Some  difficulty  does  arise  with  limitations  :n 
sizes  of  available  TPX  bulk  stock.  ExtruCeO  olate 
stock  is  onlv  availale  in  thicknesses  un  to  one 
iren  widths  un  to  13  inches)  and  extruceo  rocs 
up  to  'our  inenes  in  diameter,  "he  nrespnt  appli¬ 
cation  ronuirod  j  lens  thickness  .if  ’  •  inenes  ang 
initial)-'  it  was  felt  ‘hat  a  strong  aor.es  ive  with 
low  optical  absorption  (or  'pflpction/  would  bp 
sufficient,  "easurpmenrs  were  -ace  of  t ■?  'rans- 
11155100  at  225  GHz  of  *?X  samples  laminated  wit.i 
various  adhesives.  The  excellent  optical  'ro'er- 
ties  of  Gorlano  Optical  .Adhesive  61  mace  it  the 
best  choice,  but  reinforcement  of  the  bond  was 
necessary  for  jnacoininq.  This  was  acconnl-m.ed 
with  dowel  rods  of  TPX.  Two  TPX  sheets  were 
machined  flat,  nolished,  and  glued.  5ol:'ning 
the  surfaces  was  an  attempt  to  minimize  interface 
reflections  tiv  making  the  adhesive  lever  -ary 
thin  *he  dowel  rods  were  cooled  to  -53  I  it 
whirn  temperature  tneir  Ji.i.-eters  „erp  s-ltler 
bv  at  least  0  003  inch.  These  were  inserted  into 
holes  m  the  laminated  *PX  sheets  jnd  a  I  ’ ...g  m 
exnand  to  a  very  tight  ‘it.  After  mis  rroreoure 
machining  of  the  lens  was  successfully  ro-nletert. 
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Eccosorb  Brewster  Angle  Measurements 

Eccosorb  (TM,  Emerson  and  Cummings)  is  a  mate¬ 
rial  which  has  been  used  to  make  very  compact 
black  body  radiators  [9].  In  these,  the  high  ab¬ 
sorption  of  the  material  is  utilized  by  viewing 
the  source  at  Brewster's  angle  and  thereby  reduc¬ 
ing  reflections,  in  the  interest  of  extending  the 
use  of  this  material  to  frequencies  above  183  GHz, 
new  measurements  of  the  reflection  properties  have 
been  made. 

Determination  of  the  Brewster  angle  for  Ecco¬ 
sorb  at  163,  254,  570.5,  and  1217  ym  were  per¬ 
formed  with  a  far-infrared  laser.  With  the  equip¬ 
ment  arrangement  shown  in  Figure  2  the  reflectance 
of  plane  polarized  radiation  from  Eccosorb  was 
measured  for  several  values  of  e  in  the  vicinity 
of  the  Brewster  angle  0g. 

The  Eccosorb  was  specially  molded  into  a  flat 
sided  disc  of  diameter  2  inches  which  could  be 
easily  interchanged  with  a  standard  1 00%  reflector. 
The  polarizer  was  oriented  to  pass  radiation 
linearly  polarized  in  the  plane  of  incidence  to 
insure  maximum  penetration  into  the  Eccosorb  at 
9g.  Results  of  the  measurements  are  given  in 
Table  2. 

Metal  Grid  Polarizing  Beam  Splitter* 

A  polarizing  beam  splitter  is  often  an  impor¬ 
tant  component  in  quasi -optical  designs  for  di- 
plexers,  duplexers,  Fabry-Perot  type  filters,  or 
Fourier  *ransform  spectrometers.  One  type  is 
made  by  wrapping  thin  wires  tightly  around  a 
frame  [10],  but  these  are  difficult  to  make  ac¬ 
curately,  are  fragile  and  costly.  An  alternate 
kind  is  made  by  photo-etching  the  grid  pattern  on 
a  thin  dielectric  grid  substrate.  The  present 
work  has  been  an  attempt  to  extend  the  state  of 
the  art  in  this  medium,  so  that  larger  apertures 
would  be  available  and  good  performance  obtained 
at  high  frequencies.  Aluminum  grids  on  Mylar 
(TM,  Dupont)  have  been  constructed  with  the  fol¬ 
lowing  parameters:  strip  and  gap  widths  5  urn, 
film  thickness  6  um,  and  apertures  up  to  100  mn 
diameter.  Preliminary  tests  with  radiometers  at 
frequencies  near  200  GHz  have  indicated  better 
than  99%  polarization  efficiency  and  less  than  1% 
loss.  Fourier  spectrometric  measurements  will  be 
made  and  compared  with  theory. 
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TABLE  X 

CALCULATED  NEAP 

REFLECTANCE  TRANSMISSION 

ABSORPTION 

com. 

Utfllta 

2. 422 

570.5 

O.lOO 

0.137 

0.4S 

lutUU 

2.403 

1217 

0.100 

0.533 

0.13 

m 

2.413 

133 

0.033 

0.137 

0.74 

vx 

2.413 

570.3 

0.033 

0.337 

0.12 

m 

2.333 

1217 

0.033 

0.333 

M3 

J  -  :  ■  ? 

TABLE  II 

1 

REFLECTANCE 

(w«> 

Irevster 

163. 0 

0.03 

80.0 

236.0 

0.02 

80.0 

570.5 

0.12 

70.0 

1217.0 

0.10 

70.0 
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